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We review recent finite opacity approaches (GLV, WW, WOGZ) to the 
computation of the induced gluon radiative energy loss and their ap- 
\ phcation to the tomographic studies of the density evolution in ultra- 

relativistic nuclear collisions. 



' 1- Introduction 



Since June 2000 measurements of Au + Au reactions at the Relativistic 
Heavy Ion Collider (RHIC) at the Brookhaven National Laboratory (BNL) 
at ^/s = 56, 130, 200 AGeV (GeV per nucleon pair) have revealed a variety 
of novel multiparticle phenomena not observed previously in e~^e~' , ep, pp 
collisions at any energy nor in nuclear collisions at lower (SPS/CERN and 
AGS/BNL) energies {^/s = 17, 5 AGeV). While the bulk global observables 
such as the rapidity dependence of hadron multiplicities and transverse 
energy scale geometrically from elementary p + p collisions, the most strik- 
ing new phenomena were discovered in rare high transverse momentum 
observableflU^. The preliminary RHIC data are discussed extensively in 
the Quark Matter 2001 and 2002 proceedingj^^. 

At RHIC collider energies the hard pQCD rate of rare high pT parton 
scattering becomes sufficiently large that jets can be used to probe the dense 
quark-gluon plasma formed in nuclear collisions. In 1982 Bjorken proposed 



1 



February 6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



2 

that elastic final state energy loss of partons may "extinguish" jets in high 
energy p + p collisionJUJ. However, the elastic energy loss of partons in a 
QCD plasma of temper ature T ~ 300 MeV turned out to be too small 
{dEei/dx < 500 MeV/fmjEnSlfor jet extinction. The data on p + p jets up 
to Tevatron energies show no sign of deviations from unquenched factorized 
pQCD. Di-jet acoplanarity was proposed as another possible manifestation 
of elastic final state inter action J ^ 5 [16] pj-iQj- current RHIC data on 
nuclear collisions at 200 AGeV, th e data showed no hint of this effect. 
In some early 

studiejnnaii 

it was suggested that induced radiative 
energy loss in nuclear collisions could be much larger than the elastic en- 
ergy loss and jet quenching should become observable at least in collisions 
of heavy nuclei. W ith the development of the Monte Carlo HIJING event 
generato r ' I I predictions for the magnitude of the suppression pat- 
tern high pt hadrons were made in a first studjl^ as shown in Fig. ^ 
The middle panel shows that up to an order of magnitude suppression 
of charged hadrons was expected in the moderate pt ~ 5 GeV range. 
The input assumption was that the gluon energy loss in the plasma was 
dE/dx — 2 GeV/fm due to induced gluon radiative energy loss. The HI- 
JING model is a two component model with semi-hard and hard pQCD jet 
above pr > Po — GeV computed via the PYTHIA code and the "soft" 
beam jet fragments computed via a hybrid LUND and Dual Parton model 
algorithm. See the original articlc^^ for a detailed discussion and references. 

A selection of recent data from RHIC that confirm strong nuclear at- 
tenuation of moderate high transverse momentum charged hadron and 7r° 
production arc shown in Fig. [21 through Fig. [3 In order to appreciate just 
how remarkable and different this experimental discovery is relative to pre- 
vious lower energy data from the SPS, the quenching pattern of t:^ measured 
at RHIC is compared in Fig. 0] to the stronsenhancement of high pT 7r° 
measured in Pb + Pb at 17 AGeV at the SP$SI. 

The enhancement of high pT hadrons at SPS is an amplified version 
of the well known Cronin effect first observed in p + A collisions. As wc 
elaborate in later sections, it is due to multiple initial state interactions and 
completely masks any possible energy loss effects at the SPS. In contrast, 
at RHIC the high energy loss overwhelms the Cronin enhancement in our 
present interpretation of the data. _____ 

As discussed for example in Refs.™ it is useful to decompose the 
nuclear geometry dependence of invariant hadron distributions produced 
inA + B—fh + Xat impact parameter b into a phenomenological "soft" 



6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



JET QUENCHING AND RADIATIVE ENERGY LOSSIN DENSE NUCLEAR MATTERS 



dN,i/d7, 



Au+Au(b = 0) 




R"(Pt) 
Au+Au(b=0) 




100 - Vs w„=200 GeV 




with shadow 
and quenching 

r no shadow 
[ no quenching 

r with shadow 
[ no quenching 
no jets 

....I.... I.. ..I.. 



3 4 



Pt 



Pt 



Fig. 1. HIJING predictionsisAof the inclusive charged hadron spectra in central Au+Au 
and p + Au collisions at y/s = 200 AGeV. The competing effects of minijet production 
(dash-dotted), gluon shadowing (dashed) (assuming that gluon shadowing is identical to 
that of quarks), and jet quenching (solid) with dE/dx = 2 GeV/fm are shown. R^^(j)t) 
is the ratio of the inclusive px spectrum of charged hadrons in A + B collisions to that of 
p + P- In contrast to Au + Au, no significant quenching is expected in p + Au (or d + Au) 
since only ^ 20% initial state shadowing and Cronin effects modify the pQCD spectrum 
at high Pt ■ 



and perturbative QCD calculable "hard" components as 

dNAB{b) Npartib) dNsoft{b) , 1 d(Jhard{b) 

where Npartib) is the number of nucleon participants and Ncou{b) = 
crf^TAsib) is the number of binary NN collisions at impact parameter 
b. The nuclear geometry of hard collisions is expressed in terms of the 
Glauber profile density per unity area TAsib) — j d^T^ Ta{v)Tb{v — b) 
where TA{r) = J dz /9^(r, z) (see Fig. 6). The hard part scales with the 
number of binary collisions cx A^^^ because the probability of high pT pro- 
cesses is small and built up from all possible independent parton scattering 
processes. The soft part scales with only Npart oc A because the proba- 
bility of low transverse momenta processes is bounded by unitarity. This 
is sometimes referred to as Glauber shadowing or saturation depending 
on the calculational frame. The "soft" part is actually the overwhelming 
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Fig. 2. Preliminary STAF0 data on charged hadron quenching between central and 
peripheral Au + Au collisions at y'i = 200 AGeV is shown for peripheral event classes. 
The data are scaled by the Glauber binary collision scaling factor for each centrality 
class. Upper gray band is the expected result from pQCD scaling from p + p. Lower 
gray bands correspond to scaling with the nucleon participant number instead. Beyond 
4 GeV, charged hadrons are suppressed by 0.2-0.4 (factor 2.5 to 5) relative to naive 
binary scaling. 



bulk of the produced hadron distribution. It is expected to reflect the col- 
lective hydrodynamic properties of the produced plasma in A + A colli- 
sions. Relativistic hydrodynamics predicts specific flow patterns, such as 
hadron mass dependent transverse elliptic flo'wIS 32 33 34^ which may pro- 
vide direct constraints on the QCD equation of state. Unfortunately, low 
transverse momentum processes are not directly computable via QCD and 
many competing phenomenological models can be adjusted to fit the data 
not only at RHIC but SPS and AGS as well. 

The great advantage of RHIC over the previous AGS and SPS explo- 
rations of nuclear collisions is that the computable high px pQCD processes 
are sufficiently abundant, and that they can be used as effective "external 
probes" of the quark-gluon plasma that is produced. High pT and heavy 
mass partons are produced first, at time St l/my fm/c, while the most 
of the partons of the plasma with temperature T form and equilibrate at 
later times ~ l/<?r ~ 0.5 fm/c. Hard jets propagate along approximate 
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Fig. 3. Preliminary PHENIjJ^ data on charged and tt" hadron quenching between cen- 
tral and peripheral Au + Au collisions at yfs = 200 AGeV is shown. The pions are 
generally more quenched than the summed charged hadrons {tt + K + p). 
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Fig. 4. The striking contrast between PHENDffldata on vr" quenching in central Au + 
Au collisions at a /s = 130 AGeV compared to the (Cronin) enhancement found in Pb+Pb 
at 17 AGeV fronJSI demonstrates that jet quenching is a new nuclear phenomenon first 
seen at RHIC. 
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Fig. 5. Preliminary = 200 AGeV PHENIjP data on tt" quenching on central and 
peripheral Au + Au collisions is compared relative to preliminary PHENIX data on p+p. 
Peripheral reactions are consistent with simple binary (~ A"*/^) scaling from p + p while 
in central collisions substantial quenching is observed. 



straight eikonal lines through the plasma until r ~ i? ~ 5 fm/c. The en- 
ergy loss and transverse momentum broadening suffered by the jet prior to 
hadronization provides the tomographic handle to probe the opacity of the 
evolving plasma. Since non-central collisions have a well known geometrical 
asymmetry, the azim uthal d istribution of final high px hadrons provides 
even more informatiorP^EHl. This is illustrated in Fig.[7|by the normalized 
average thickness of nuclear matter as a function of the angle with respect 
to the reaction plane that a fast parton sees on its way out. The azimuthal 
distribution can be Fourier decomposed into 



dN{b) _ 1 dn{b) 



+ V2 (y, b, pt) cos(2(0 - ^fc)) 



(1 + 2wi [y, b, Pt) cos{(f> - (/){,) 
••) , 



(2) 



where (jjb is the azimuthal angle of the impact parameter b. In Fig. |S1 the 
large elliptic flow component V2 is seen to extend up to the highest pT 
in non-central collisions at mid-rapidity. Another critical new experimen- 
tal test for jettiness at moderate pt ^ 5 GeV is provided by two parti- 
cle correlation. The data from both STAR and PHENIX show clear evi- 
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Fig. 6. Illustration of key aspects of the relation between the geometry of nuclear colli- 
sions and the participant Npart{b) = 2 J d^rTA(b- r){l - exp(-o-PPTA(6))) ~ < 2A 
and collision N^allib) = o-PPTAA{b) ~ j4*/^ nun rber at a fix ed impact parameter. The 
observables dN'^^ /dr], dET/dr], V2(pt) (see RefsEZEHESEOI) used to constrain the ge- 
ometry experimentally are also illustrated. 



dence for back-to-back correlations predicted by pQCD (and checked via 
the PYTHIA generator). However, there appears to be a difference between 
the correlation pattern of charged and neutral hadrons as seen in Fio-s IHl 
and El The STAR data are in accord with predictions in Refs!' 
that the two jet structure should be quenched by an order of magnitude 
relative to jet correlations 'm p -\- p due to the high opacity of the plasma. 
Similar jet structure in high pTneutral hadron correlations is also found in 
the preliminary PHENIX datcP. Whether the suppression of away-side jet 
seen in STAR data also manifests in the PHENIX experiment has yet to 
be addressed. 



6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



1.4 

1.2 

1.0 

^ 0.8 
n 

a 0.6 
0.4 
0.2 
0.0 



Au+Au 



Woods-Saxon geom. 
Sharp cylinder geom. 



b = 2fm 




0.0 



2.0 4.0 
^ [rad] 



6.0 



Fig. 7. The average normalized optical depth (R from Fig. 2 of RefEHl) seen by a parton 
propagating through a plasma in azimuthal direction (f> for different impact parameters 
and different density profiles. 
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Fig. 8. Preliminary = 200 AGeV STAR on high pT elliptic flow is shown. A value 
of V2 = 0.2 corresponds to a 2-to-l azimuthal asymmetry of hadrons. 
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Fig. 9. STAR data on high azimuthal correlations of charged hadrons at 200 AGeV 
from R,cf>^ showing clear evidence near side and away side correlations characteristic of 
two jet production in elementary p + p. However, in Au + Au reactions the away side 
jet correlations are reduced in central collisions. The Back-to-Back jet correlations are 
attenuated by an order of rnagnit u de in the most central collisions as expected if the 
plasma produced is opaqu o I I I 

For our review, the most important point about both data sets is that 
narrow azimuthal near side and away side correlations do exist. This is a 
necessary (though not sufficient) conditions to enable us to interpret mod- 
erate pt hadronic quenching patterns as due to jet energy loss. 

In this review we focus on the rare hard pQCD tail of the hadron yields 
where all these interesting new nuclear phenomena have been discovered 
at RHIC. The central element of the theory that will be elaborated in the 
following sections is the predicted dependence of the induced QCD radia- 
tive energy loss on the jet energy, the plasma density, and the expansion 
properties of the plasma. We restrict the discussion to two complementary 
approaches, the GLV and WW/WOGZ formalisms, that are best suited in 
our opinion to applications to nuclear collision problems where the opacity 
is necessarily finite and the energy range of minijets accessible experimen- 
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Fig. 10. Preliminary PHENIX datS on high px azimuthaEHl correlations of neutral 
pions also show clear jet structure amid an elliptic flow component. 



tally is typically 3—15 GeV. 

In the absence of final state interactions we recall the well-known lowest 
order invariant pQCD differential cross section for inclusive p + p ^ h + X 
given by 



(pp 



= K 



abed 



dXadXbfa/piXa, Qa)fb/p{Xb, Qb) 



X — ao - 

dt 



cd)- 



D 



(3) 



where Xa = Pa/ Pa, Xb = Pb/Ps are the initial momentum fractions carried 
by the interacting partons, Zc = Ph/Pc is the momentum fraction carried 
by the final observable hadron, ,fa/p{xa,Qa) are the parton distribution 
functions (PDFs), and Dh^c{zc, Q^) the fragmentation function (FFs) for 
the parton of flavor c into h. 

The phcnomcnological K factor is introduced to mimic next-to-leading 
order (NLO) corrections. One finds that Eq.© tends to over-predict the 
curvature of the inclusive hadron spectra in the pt < 4 GeV range. This can 
be partially corrected via the intrinsic fc^-smearing of partons associated 
with vacuum radiation and generalized parton distributions fa{x, kr, Q^)- 
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For the corresponding modification of the kinematics in J^Jl in addition to 
the / (Pk^ J cPk^ (■ ■ ■ ) see Refs.^^ The generalized parton distributions 
are often approximated as 

U{x, kr, Q^) « Q^)g{kT), gikr) ^ ^^^^^ , (4) 

where the width {k"^) of the Gaussian is related to initial state vacuum 
radi ation. Discussion of some qualitative features of (fc|^)pp is given in 
Refs.^^^For comparison to experimental data the reader is referred to 
Fig. ^] shows that this parton model approach provides a 
good description of the high pT data on hadron production in the range 
above a few GeV at all center of mass energies. 

In order to calculate the effects of parton energy loss on the attenua- 
tion pattern of high pT partons in nucle ar collisions, we must modify the 
free space fragmentation functionJ^^^21 Energy loss of the parton prior 
to hadronization changes the kinematic variables of the effective fragmen- 
tation function. In the first approximation, this effect can be taken into 
account by replacing the va cuum fragmentation functions in Eq.© by ef- 
fective quenched one^SESl 



I _ Ph ^ Ph 

Pc~ AE,{p,, 0) ' ' AE,{p,, ct>)/Ng ' 



(5) 



where z^, Zg are the rescaled momentum fractions. The first term is the 
fragmentation function is of the jet c after losing energy AEc{pc, 0) due to 
medium induced gluon radiation. The second term is the feedback due to 
the fragmentation of the Ng{pc, (f)) radiated gluons. The modified fragmen- 
tation function satisfies the light-cone sum rule j dzc ZcD'i^^^{zc,Q^) = 1. 
Eq.(|SJl takes into account the dependence of the energy loss on the par- 
ent parton energy and the possibly azimuthally asymmetric region of high 
density nuclear matter. 

In the above first approximation, only the average val ue of the energy 
loss is used. One can also include multigluon fiuctuationJ^^EZI of the en- 
ergy loss via an energy loss distribution P(e, E) where e = coi/E is the 
fractional energy loss of a jet of energy E in the rest frame of the plasma. 
The mean energy loss in the first approximation is related to P via 

de P{e,E)e = AE/E . (6) 
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Fig. 11. Comparison of the pQCD modeBHwith intrinsic (fc|,> = 1.8 GeV^ and GRV 98 
PDFs to p+p data at y/s = 540, 1800 GeV. The normaUzation K factor simulating higher 
order corrections is fit at each energy. A ±25% error band is also included. 



The invariant hadron distribution attenuated by fluctuating energy loss 
m A + A collision is then given by 



AA 



(Pp 



TAA{b) ^ j dxidx2 J d^Kd^kb gika)g{kb) 

abed 

X SA{Xa,Ql)SB{xb,Ql) fa/A{Xa,Ql)fb/B{xb,Ql) 



(7) 







February 6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



JET QUENCHING AND RADIATIVE ENERGY LOSSIN DENSE NUCLEAR MATTER13 

where z* = 2:^/(1 — e). Here, T4yi(&) is the Glauber binary collision profile 
density at impact parameter b. It is found numerically , that fluctuations 
tend to reduce the magnitude of the attenuation by a factor about two at 
RHIC. It is important to include this effect when inverting the attenuation 
pattern in jet tomography to determine the initial plasma density. In B + A 
reactions isospin effects can be accounted for a nucleus with Z protons and 

N neutrons by fa/A{Xa,Ql) ^ (Z/A) fa/p{Xa,Ql) + (N/A) fa/n{Xa,Ql)- 

Nuclear modifications to the PDFs are estimated in our applications by 
using the shadowing function SA{xa,Qa) proposed by EKS'9&^^. 

2. GLV Theory of Radiative Energy Loss 

In this section we turn to the theoretical problem of computing the energy 
loss of a fast parton penetrating a finite, expanding quark-gluon plasma. 
High pt many body pQCD is a new frontier of research at RHIC and even- 
tually LHC. The needed new theoretical development is the non-Abelian 
analogue of the radiative energy loss theory familiar from classical E&M. 
An important practical problem in QCD is that there exists no external 
beam of isolated colored quarks or gluons, and that the jets are always pro- 
duced in hard processes just before the plasma is created. In addition, the 
nuclear medium has a small dimension compared to the jet fragmentation 
coherence length. T herefore, t he basic formation time physics of Landau- 
Pomeranchuk-Mig (LPM) is expected to lead to strong destruc- 

tive interference effects that have to be carefully tak en into acco unt. 

In this section we review one of those approache d I I I I developed 
by Gyulassy-Levai-Vitev (GLV) based on an algebraic reaction operator 
formalism. The expansion parameter is the opacity x — L/ X = apL of 
the system and the result is presented to all orders in powers of x or 
equivalently to all twist parton-parton correlatio ns. Othe r approac hes rely- 
ing on asympto tic techniques include BDMS^^ 5 7 | 58|59| 26O 61 62^31^ 
SwEESESEZI have been reviewed elsewher^SSl in the second part of this 
report the twist expansion approac d^^ l ^^ l ^ ^ for parton energy loss in nu- 
clear matter is reviewed. 

2.1. The GW Plasma Model 

The GLV approach is built around a simple model of multiple scatter ing 
in a plasma formulated in gW22J However, the results of GL\E2ES1 are 
more general since even out of thermal equilibrium the effective in-medium 
interactions are of finite range R ~ M^^- Consider the sequential elastic 



February 6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



14 

scattering of a high energy (jet) parton in the random color field produced 
by an ensemble of m static partons located at Xj = (zi,x_Lj) such that 
Zi+\ > Zi and {zi^i — Zi) » where ^ is the color screening mass in 
the medium. As a simplified model of multiple scattering in a color neutral 
quark-gluon plasma, we assume a static Debye screened potential for each 
target parton: 

V^i"(q)=ff(i;")c,c'-2^e-"i- , (8) 

where T" is a dj-dimensional generator of SU {N) corresponding to the rep- 
resentation of the target parton at x^. The initial and final color indices, 
c, c', which refer to the target parton, are averaged and smnmed over when 
computing the ensemble averaged cross sections. With V-^ cx the ensem- 
ble averaged potential vanishes everywhere, {V") oc TrT" = 0. However, 
since 

TrT^T^ = 5ab5ii{di/dA)C2i , (9) 

the diagonal mean square fluctuations and the cross sections are finite. 
Recall that for SU{N) the second order Casimir, C2i = {N'^ - 1)/2N = Cf 
for quarks in the fundamental {di = N) representation, while C2i = N = 
Ca for gluons in the adjoint (dj = N'^ — 1 = (Ia) representation. 

In this potential, each scattering leads on the average to only a relatively 
small momentum transfer — {<^i■,<lzi■,^.l.i) with each component being 
much less than the incident energy, Eq . The assumption that the potentials 
are static is approximately valid in a high temperature plasma of massless 
quarks and gluons in the following sense: As T ^ oo, the effective coupling 
g ^ (albeit very slowly). The perturbativc Debye screening mass ji ^ 
gT limits q±^gT. The typical thermal energy Et ~ 3T of the plasma 
constituents is therefore large compared to /U. Consequently, the average 
energy loss per elastic collision, —q^ « —q^ « q\/2ET oc g^T, is g times 
smaller than the average transverse momentum transfer. 

Because we are interested in relatively low momentum transfer scatter- 
ing [Kqcd <^ q± ^ gT <^ T), the spin of the partons can be neglected. 
The jet parton is allowed, however, to be in an arbitrary d-dimensional 
representation of SU{N) with generators, T'^, satisfying T^T" = 6*21^. 

The Born (color matrix) amplitude to scatter from an incident four 
momentum pf_-^ to in the potential centered at Xj is then given by 



Mi{pi,pi-i) = 27r^(p0 -pti)A,(q,)e-^' 



(10) 
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where qj = — Pi_i, and Ai is shorthand for 

A;(qO = r'^A"(qO = -2igEoT'^Vt{q.) . (11) 

The differential cross section averaged over initial and summed over final 
colors of both projectile and target partons reduces to the familiar form for 
low transverse momentum transfers: 

da,/dq^, ~ C*,;-^ — ■ — — , (12) 

where the color factor is 

C,^^Tr{T-T'')Tr{T^T^)^C2C2^/dA • (13) 

For SU{i), the number 2Ci gives the usual color factors 4/9,1,9/4 for 
qq, qg, gg scattering respectively. In our notation, the angular distribution 
is given by 

da,/da, = ^Tr|^,(q,)|V(47r)2 . (14) 
2.2. GLV Formalism 

In Refsl^Ella systematic recursive graphical technique was developed and 
translated into an algebraic operator method. The goal was to compute 
medium induced gluon radiation amplitudes of the type shown in Fig. 1121 
The exponential growth of the number of graphs with the number of in- 
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Fig. 12. Induced radiation amplitude^sJ contributing to fif th o rder and higher orders in 
the opacity expansion of QCD energy loss in the GW modeC2l The crosses denote color 
screened Yukawa interactions on a scale ^. The blob is the initial hard jet amplitude. 
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teractions makes it very tedious to go beyond order three in opacity x or 
twis t § (2 + 3x2) from the parton-parton correlations in the medium. In 
GLV ^^ I Otf l the combinatorial problem is solved by summing the inclusive 
radiative gluon distributions recursively. The first step in the approach is 
to compute the three direct (single Born) and four surviving virtual (con- 
tact double Born) diagrams shown in Figs . II 31 and 1 1 41 that contribute to the 
first order in opacity induc ed r adiation. Detailed instructive calculation of 
these amplitudes is given iJ^. 
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Fig. 13. Three first order (singe Born) direct amplitudes that s erve to define the b„ 
component of the reaction operator in Eq. I21i are derived in Refsinimi 



For scattering of a jet off of n scattering centers located at depths Zi in a 
homogeneous medium of large area fJ,'^A± 3> 1, we can write the inclusive ra- 
diative gluon "probability" , P„ (k, c) , as a sum over products of partial sums 
of amplitudes and complementary complex conjugate amplitudes'^ Ev- 
ery term in the sum contributes to the same 0(5^"^^) in the strong coupling 
constant g, where comes from the radiation vertex and (g*)" is related to 
the n elastic scattering cross sections. The average value of n is th e opacity, 
X — {n) — L/X, of the medium. It turns out surprisingljE^I^ that A is 
the mean free path of radiated gluons rather than the jet itself, as we shall 
show. The partial sums of diagrams at order n in the opacity expansion can 
be conveniently expressed in a tensor notation and constructed by repeated 
operations of 1, _D,j, or Vi corresponding to no, direct, or virtual interactions 
at scattering center i: 

n 

^ii...i„(a;,k, c) = Y\_ {^Q.im + ^i,i^Drn + S2,i„,Vrn^ Go{x,'k,c) . (15) 

m—1 

Here Go is the initial hard {q,g} + g color matrix amplitude, x = /p'^ is 
the gluon lightcone momentum fraction, k is the gluon transverse momen- 
tum, and c is its color index. 

In the inclusive probability each class contracts with a unique comple- 
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Fig. 14. Diagrams with two momentum exchanges that contribute to first and second 
orders in opacity. In the contact (z2 = zi) limit the (double Born) virtual (contact) 

Note 

of the radiation 



amplitudes define the Vn components of the reaction operator in Eq. I21L 
that two of the diagrams vanish due to zero measure J^^ dz 

interval and two are topologically indistinci^^ESl 
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Elastic 



(Zl) 



(Z3) 




Inelastic 




Fig. 15. Example of graphs constructed via Di,Vi that contribute to the 4-th order 
in opacity (or twist 2+8 in-medium partoii-parton correlations) in elastic and inclusive 
inelastic final state interactions. The longitudinal depths of active scattering centers are 
denoted by Zi and inactive (created with 1;) by (zi). The form of Di,Vi depend on the 
process type but the tensorial bookkeeping of partial sums of amplitudes is the same. 



mentary class, 

P„(x,k) = A'i-'"(,x-,k,c)Ai....„(a:,k,c). (16) 

The complementary class is constructed to satisfy the fixed opacity power 
oc x" requirement as 

n 

A'^---'-{x,k,c) = Gl{x,Kc)]l(^So,i^Vi + 6,,i^Dl + 62,i^) • (17) 

m=l 
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Fig. 1151 shows an example of how this formahsm works at 4th order in 
opacity for elastic and inelastic inclusive distributions. 

Detailed diagrammatic calculations that illustrate the color and kine- 
matic modifications to the amplitudes for both direct and double Born mo- 
mentum and color exchanges with the propagating jet-|-gluon system are 
given m Refs.™ Direct single momentum transfer interactions enlarge 
rank n — 1 class elements as follows: 

-Dn.4ij ...i^ j (x, k, c) — (a,j + S*,! + i?n).4ij (x, k, c) 

= tt^^wAj^ (x, k, c) -j- e ^ ^ ".Aj^ (x, k q^,[c, a^]) 

^ ^ A'„(Ai...,„_i) 



B„e'"«^"[c,a„]rei(Ai...»„_J , (18) 

where B„ = H — C„ — k/k-^ — (k — q„)/(k — q„)^) is the so-called Bertsch- 
Gunion amplitude for producing a gluon with transverse momentum k in an 
isolated single collision with scattering center n. The momentum transfer 
to the jet is q„. The notation a;„ = (k — q„)^/2w is for a gluon with energy 
uj and a„ is the color matrix in the dji dimensional representation of the 
jet with color Casimir Cr. Ny — X]m=i counts the number of virtual 
interactions in Ai-^...i^_-^. Tei{Ai^...i^_-^) is the elastic color factor associated 
with all n — 1 momentum transfers from the medium to the jet line. 

Unitarity (virtual forward scattering) corrections to the direct processes 
involve the sum of four surviving double Born contact diagrams in Fig. 1131 
that enlarge rank n — 1 classes via the action of 

Vn = ^\i'^A + Cr) - anSn - a„-B„ = -a„z)„ - i(CA - Cr) . (19) 

This key op erator r elationship between direct and virtual insertions discov- 
ered in 

RefsPSa 

reduces the problem to recursive algebra. More specifi- 
cally: 

KAi- i„_i(a:^,k, c) = ^-^(Ca + Cr) - a„S'„ - a„i?„^ Ai - j„_i (a;, k, c) 

i \ ^20) 

The tensorial bookkeeping of classes of diagrams makes it possible to 
construct the distribution of radiated gluons in the case of n interactions, 
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P„(x,k), recursively from lower rank (opacity) classes via a "reaction" op- 
erator 

p„ = A''■■■'-'RnA^,...^^_, , Rn = dId„ + K + V,] . (21) 
Using the key identity 119(1 . the reaction matrix simplifies to 

Rn = (Dn - an)\bn ~ «„) " C (^n + B„)t(^„ + B„) - Ca ■ 

The next major simplification occurs because both S and B involve the 
same gluon color rotation through ijca„d_ xhis fact reduces the color alge- 
bra to simple multiplicative Casimir factors in the adjoint representation 

A'^-'-'{SlSn " Ca)A,. = Ca (F„-i(k- q„) - i^„_i(k)) 

= C^(e''i"-^-l)p„_i(k) , (22) 

where b = jV k is the transverse momentum shift operator. We have proved 
in Refs.E2E3that 

A''-'-'BiBnA,...^„_,=0 . (23) 

The interference term is also found to be color trivial (containing only pow- 
ers of the Casimir in the adjoint representation Ca) and can be expressed 
in recursive form: 

2-ReA''-'-'BlSn^,,...,^_, = -2CaB„- (Rc e-^""^"e^i"-^I„_i) . (24) 

Eas. H22|l - H24|l present a clear proof that only the gluon mean free path 
which is proportional to Ca (or equivalently the gluon elastic scattering 
cross section) appears in the induced non-Abelian bremsstrahlung formulas. 
The transverse vector amplitude I„ above obeys a recursion relation 

I„ = CAfe'("«-"")^"e^'»"-'^-l)l„-i-<5„,iC^CflBie'"°^^ , (25) 



where Iq = — CflHe*""^", and the Gunion-Bertsch amplitude can be written 
as 

g^gic^o^i ^ _ |^gi(^o-^i)2igiqi-b He^^o^i ^ (26) 

For n > 1 Ea. H25|l can be solved in closed form 



n — '^R'^A 



m—1 



' (27) 



where the product is understood as ordered for left to right in decreasing 
order in operators labeled by m. 
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The inclusive radiation "probability" is then found to obey the soluble 
recursion relation 

P„(k) = CA(P„-i(k - q„) - P„_i(k)) 

- 2CaB„ • (Re e-^'^-^"e^'»"-^I„_i) + 5„,iCaC^j|Bi|2 . (28) 

The initial condition for this recursion relation is the initial hard vertex 
radiation amplitude without final state interactions that is given by Pq = 

The solution to the problem for any order n in mutiparticle interactions 
can therefore be expressed in closed form as 



P„(k) = -2CRClReJ2 { n (^'"''^ - 1) \ ■ ^''''■^e" 



X \ n (e*('^°-"''")^'"e''»"'-^ - 1) \ 0H(e' 



(29) 



This expression is very general and is suitable for any distribution of interac- 
tion centers as well as any distribution of Zi dependent transverse momenta 
exchanges that can arise due to expansion of the medium. Therefore, it 
can be directly applied to realistic expanding plasmas where the distance 
between adjacent centers (the mean free path) varies along the medium as 
does the local screening scale, iJ.{z). This form is also well suited for possible 
future Monte Carlo implementation for arbitrary q_i,Zi medium ensemble 
averages. 

The final complete solution to the inclusive induced gluon radiation 
valid to all orders in opacity can be expressed in terms of the following 
infinite series: 



_dNs 
' dx 



(Pag{zi) 



- ag{zi)5'^{(ii) 



Pnizi, - ■ ■ , Zn) I —2 C(i_.,. ^n) ' ^ '^(m+1,--- ,n)(m,--- ,n) 



cos 



^'^{k,-,n)^Zk \ - COS I ^W(fc,...,„)A^fe 
, fe=2 / \fe=l J 



(30) 
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where ^2 = is understood. The notation is defined as follows: 

TT^Jl P (k-Qn -q,, q»„) 

k^' ^(n..-.™) (k-q,^-q,^ q,^)2' 



_ (k - q„ q„)" 

-^(m,--- ,n) ~ 



(31) 



2xE 

The infinite opacity series can be understood as a sum over all twist parton- 
parton correlations, (• • • {2nFF) ■ ■ ■), in the nuclear matter. We emphasize 
that Ea. (|3()(l is not restricted to uncorrelated geometries as in Refs.l'^^^H^. 
It also allows the inclusion of finite kinematic boundaries on the q^ as well as 
different functional forms and elastic cross sections (Jg{i) along the eikonal 
path. The fi rst and second orders in opacity (twists 2 + 2, 2 + 4) have also 
been checkecJ^^^ through explicit calculations of the corresponding direct 
and virtual cut diagrams. The ~ k divergences, naively present in 

the propagators C(^,^,^...i^^),B(^,^,^...i^fQ^j^...,^^) in Eq.ljSOI), are canceled by 
the interferences phases (^(m. ■■■.«)■ Ejq.lljnjl also provides further insight in 
the LPM effect in QCD at a microscopic/diagrammatic level by keeping 
track of the destructive interference among all pairs of amplitudes to all 
orders in opacity. 

For particular models of the target geometry one can proceed further 
analytically. For a homogeneous rectangular geometry, the average over the 
longitudinal target profile with 

Pnizi,- ■ ■ ,Zn) ^ n\ Pq 0{L - Zn)9{Zn - Zn-l) ■ ■ ■0{Z2 - Zi) , (32) 

where the mean density is po = Ns/LA_i_, leads to an oscillatory pattern 
that is an artifact of the assumed sharp edges. Here A± is the transverse 
area of the interaction region and is the total number of scattering cen- 
ters. A somewhat more realistic model may be an exponential longitudinal 
distribution of scattering center separations 

P„(.i,---,..) = n^^^^e-^^^/-=(") . (33) 
Le(n)A± 

This converts the oscillating formation phase factors in Ea. H3Uf) into 
Lorentzian factors 

f / " \ iV" ™ 1 
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In order to fix Lf^{n), we can impose (zj, — zq) — kL/{n + 1). This 
constrains^^ ^3 54 ^ /./(n + 1). 

2.3. First Order Non-Abelian Energy Loss 

The simplest and fortunately dominant (a posteriori ) applicatioiJS2ESl of 
our general solution to the energy loss problem was to calculate numerically 
the total radiated energy loss as a function of jet energy E, plasma depth 
L, and the typical transverse momentum transfer fi. In the absence of a 
medium, the gluon radiation associated with the parton jet (in the small 
X approximation, where we do not distinguish between quark and gluon 
parents) is distributed as 

where x = /E^ w (^/E, and Cr is the Casimir of the jet in the (Ir di- 
mensional color representation. The differential energy distribution outside 
a cone defined by k'^ > fjl^ is given by 

^.^i^log^, (36) 
ax TT ^ 

where the upper kinematic limit is k^ax — miii[AE'^x'^,AE'^x{l — x)] . 
The energy loss outside the cone in the vacuum is then given by 

A£;(«) = ?^£;iog-. (37) 

TT /i 

Beyond the small x approximation, the splitting functions Pgq{x) and 
Pgg{x) can give rise to some small corrections. While Ea. (j37(l overestimates 
the radiative energy loss in the vacuum (self-quenching), it is important 
to note that Ai?'^°) /E ^ 50% is typically much larger than the medium 
induced energy loss. However, the vacuum energy loss is included in the 
DGLAP evolution of the fragmentation functions Dh/c{z, Q^). 

To compute the medium induced radiation we focus on the density of 
the Ns scattering centers given by Ea. H33|l . Averaging over the momentum 
transfer qi via the color Yukawa potential leads to a very simple first order 
opacity result for the x ^ 1 gluon double differential distribution 



^ dxdk^ ^ dxdii^ Xg Jq ^ 7r(q^ + fi'^y 

2k.q,(k_q,)^X2 (38) 



16x^E^ + (k- qi)4L2 
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where the opacity factor L/\g = Nga^,^^ /Aj_ arises from the sum over the 
Ng distinct targets. Note that the radiated gluon mean free path Xg = 
{CA/CftjX appears rather than the jet mean free path. The upper kinematic 
bound on the momentum transfer is q^ax = s/4 ~ S-E/i and 1//^^^^ = 
1//Lt^ — l/in^ + (Zmax)- For SPS and RHIC energies this finite hmit cannot 
be ignored. Numerical results comparing the first three orders in o pacity 
corrections to the hard distribution Eq.JSHJ were presented in Refs.'"^'^ 
and the opacity series was shown to converge fast is dominated by its first 
order term for realistic nuclear opacities. Fig . II 61 illustrates the convergence 
properties of Eg . 1)30(1 in the example of the mean energy loss AE, as well 
as the quadratic dependence of AE on the size of the nuclear matter for 
static media. 
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Fig. 16. The radiative energy loss of a quark jet with energy -Bjet = 5, 50, 500 GeV (at 



SPS, RHIC, LHC) is plotted as a function of the opacity L/\g {Xg 



1 fm. 



: 0.5 GeV) 



for a static medium from RefE2E3l Solid curves show the first order in opacity results. 
The dashed curves show results up to second order in opacity, and two third order results 
are shown by solid triangles for SPS energies. 



In Ref. we considered various asymptotic limits of the first order en- 
ergy loss in 1+lD and 1+3D expanding plasmas. In the case that the mean 
density decreases as 



p{z 



Po 



(t 



(39) 
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due to a longitudinal Bjorken expansion, it is possible to obtain a closed 
analytic formula under the strong assumption that no kinematic bounds 
need to be considered. For a hard jet penetrating the quark-gluon plasma, 

^ ---/(Z(x,r)), (40) 



dx TT J^^ A(r 

where a; ~ cj / E is the momentum fraction of the radiated gluon and the 
formation physics function f{Z{x,T)) is defined^ to be 

/•oo 7 

fix, r) = / — ^ [1 - cos ( uZ{x, r))] . (41) 
7o u{l + u) 

With Z{x,t) = (r — To)fi^{T) /2xE as the local formation physics param- 
eter, two simple analytic limits of Ea. (|41|l can be obtained. For x ^ Xc — 
/i(ro)Vo=' L^- — /2E = LiP{L)/2E, in which case the formation length is 
large compared to the size of the medium, the small Z(x,t) limit applies, 
leading to f{Z) « 7rZ/2. The interference pattern along the gluon path 
becomes important and accounts for the the non-trivial dependence of the 
energy loss on L. When x <C Xc, i.e. the formation length is small compared 
to the plasma thickness, one gets f{Z) « logZ. In the x ^ Xc limit-^*'-, the 
radiative spectrum Ea. H4U|) becomes 



dx ^ 2(2 - a) A(ro) x 
In the opposite x <^ Xc limit we have 



(42) 



dx ^ ^(3 - a) A(to) ^ 2xE ^ ' 

and the intensity spectrum has an integrable divergence at a; = 0. Higher 
orders in opacity do change the x-dependence of the radiative spectrum 
relative to Eqs.(g2l ^^'^ ESIj but this change is smah^^ ^for X i?^ ~ 
5 fm. 

The mean energy loss (to first order in x) integrates to 

The l ogarit hmic enhancement with energy comes from the Xc < x < 1 
regiorP^li^. In the case of sufficiently large jet energies (E — > oo) this term 
dominates. For parton energies < 20 GeV, however, corrections to this 
leading log 1/xc expression that can be exactly evaluated numerically from 
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the GLV expression and are found to be comparable in size. The effective 
llE / E in this energy range is approximately constant. 

In order to study non-central collisions and azimuthal anisotropy we also 
solved the analytically tractable case of a sharp expanding elliptic cylinder. 
We approximate the assumed 0o independent screening /x(t) « gT{T) = 
2(/9(t)/2)1/3 since 5 ~ 2 and p = (16C(3)/7t2)T3 ~ 2T^ for a gluon plasma. 
We define r(0o) as the escape time for the jet to reach the expanding elliptic 
surface from an initial point Xq = (xo,?;o) in the azimuthal direction 

{xq +TCOS(j)oY {yo + Tsin^o)^ ^ 
(R. + v^tY + {Ry + VyT^ ~ ■ ^ > 

We take a;((/)o) = 2t((/>o) (p(r(0o))/2)^^"^ to estimate an upper bound on 
the logarithmic enhancement factor. Performing the remaining integrals 
one gets: 

A..m/.x QCRaldN9 E f°°^ 1 1 



xeii- 



dy uji4>a) Jo Rx + v^t Ry + VyT 
(xo + TCOS(/>o)^ , (2/0 + Tsin^o)^ 



{Rx + V^tY {Ry + VyTY 

^ 9 C^dN^ logT+;i;7y E ^ 

4 RxRy dy ttx - ay io{<l>o) 

where = v^/Rx, ay = Vy/Ry. This expression provides a simple analytic 
generalization that interpolates between a pure Bjorken I+ID expansion 
for small a^.yT, and a 3+lD expansion at large ax.yT. 

In the special case of a pure Bjorken (longitudinal) expansion with Vx = 

Vy = 0, 

In this case, the energy loss depends linearly on t((^). 



2.4. Centrality and Rapidity Dependence of the Cronin 
Effect at RHIC 

In the subsections that follow we will describe some of the important appli- 
cations/predictions of the GLV reaction operator formalism to initial and 
final state multiparton interactions. We emphasize that these calculations 
can provide complementary information about the properties of cold and 
hot nuclear matter such as the transport coefficients [J? jXq-, \J? I\ and the 
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effective initial gluon number density and energy density created in rela- 
tivistic he avy ion reactions. 

In Refs! '^ 1 * I the explicit solution for the transverse momentum distri- 
bution of partons that have traversed cold nuclear matter has been found 
using the GLV approach: 



II ^ ri=0 



1 daei{R,T) 



■^"^ (g)e"'^'"^''ii - 1 



d^N^{kn,k^ 
c?fc||C?^k^ 



(48) 



This leads to nuclear induced broadening (Afc^)pA = /^^xC for opacity 
X ~ (L) /A and typical transverse momentum squared /z^. Beyond the naive 
Gaussian approximation harder fluctuations along the projectile path lead 
to a logarithmic enhancement ( Ak'^ )pA, i-c f = ln(l-|-cp|n). The question of 
the effective longitudin al mom entum shift associated with this broadening 
has also been addressecP^^l^ 

1 



-Afc|| - — . (49) 



To implement initial state elastic and radiative energy loss we focus on the 
large ~ partonic subprocess a + 6 — > c + d, where ka, kb are the initial 
momenta involved in the hard part da/dt of Eq.®. If partons a and b have 
lost fractions Cq, { a = a,b) of their longitudinal momenta according to a 
probability distributions Pa(e), then ka — ka/l — Sa and 



(50) 



at asymptotic t = — oo . Ea. H50() provides a simple modification to the 
factorized pQCD hadron production formalism. For bremsstrahlung pro- 
cesses, Paii) are sensitive to multiple gluon emissio For the simpler case 
where one considers only the mean energy loss, P{e) — 6{e — (Afco)/fco). 
More specifically, for the elastic longitudinal shift that we consider here, 
Pc.ie)=5(e-n^Xc.^/2kl.) and 



II. 

fa/piXa: Q^) fa/p ( + ^ -, ) 6'(ia < 1) . (51) 



The observable effects of Eas. H5()ll and H51|l can be very different for valence 
quarks, sea quarks, and gluons due to the different x-dependence of the 
PDFs. 
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Fig. 17. Rapidity dependence of the Cronin effect in d + Au reac tion s at \/^jvjv ~ 
200 GeV with and without antishadowing or EMC effect from RefBSl The result of 
switching off elastic energy loss is also shown via the upper bands for y = ±3. 



A natural application of these results is in the quantitative description of 
the Cronin effect. Perturbative calculations that include initial state parton 
broadening remain its only successful explanation to date'^^ 4 7o 7 6 77 78 
Comparison to low energy data on proton-nucleus collision j^^l^'^' allows 
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US to extract the transport coefRcients of cold nuclear matter, fJ^ /\q = 
0.06 GeVVfm and fx^/Xg = 0.14 GeVVfm. These estimates include the 
effective elastic initial state energy loss and are ~ 25% bigger rela tive to 
the case when the longitudinal momentum shifts have been neg lectejm. 
The transport coefRcients of cold nuclear matter change only slightljE^ 
with ^/s and can be used to predict the centrality and rapidity dependence 
of the Cronin effect in d + Au reactions at RHIC. 

Fig. II 71 shows the computed rapidity dependence of the Cronin effect in 
d + Au reactions at RHIC. At midrapidity, y = 0, in the computed trans- 
verse momentum range, we find slight Cronin enhancement RdAu max = 1-2 
relative to the binary collision scaled p+p result even with initial state elas- 
tic energy loss. Including strong antishadowing (shown with dashed lines) 
leads to RdAu max = 1-3. The effect of antishadowing becomes increasingly 
important at higher px- Experiments at RHIC can thus help constrain the 
poorly known nuclear modification to the PDFs for gluons. At forward 
y = +3 rapidity (in the direction of the deuteron beam) the Cronin inter- 
cept {RdAu = 1) and maximum RdAu max are both shifted to slightly higher 
Pt- The most distinct prediction at forward rapidities is the significantly 
larger Cronin enhancement region extending to highp^. This is understood 
in terms of the softening of the hadron spectra away from midrapidity, as 
predicted by perturbative QCD. Steeper spectra tend to enhance the effect 
of otherwise identical transverse momentum kicks. At backward y — ~3 
rapidity there is no enhancement since the nucleus does not scatter mul- 
tiply on the de uter on. This region is shown to be sensitive to the initial 
state energy losM For further details on the systematic imderstanding of 
the Cronin effect versus ^/s the reader is referred to Refs.^^^^^. Recently it 
has been suggested by models based on gluon saturation, final state hadron 
absorption, and coherent nucleon scattering, that hadron cross sections in 
d + Au reactions will be suppressed by about 30% for charged hadrons rel- 
ative to the binary collision scaled p -I- p result. Experimental data at RHIC 
will make possible a critical test of the validity of different theoretical mod- 
els. 



2.5. Jet Tomography of d -\- Au and Au + Au at SPS, 
RHIC, and LHC 

Some of the most important application^^l^ of the GLV results on 
medium induced gluon bremsstrahlung are related to jet tomography, the 
study of the properties of matter through the attenuation pattern of fast 
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particles that propagate and lose energy as a result of multiple elastic and 
inelastic scatterings. Under the assumption of local thermal equilibrium all 
relevant scales in the problem, such as the Debye screening scale fi, can be 
related to dN^ /dy - the initial effective gluon rapidity density of the bulk 
soft background matter (the quark-gluon plasma). 
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Fig. 18. The suppression/enhancement ratio fi a a (vt ) {A = B = Au) for neutral pions 
at Vijvjv = 200, 5500 GeV from Refi^ Solid (dashed) lines correspond to the 
smaller (larger) effective initial gluon rapidity densities at given y/s that drive parton 
energ y lo ss. Data on tt" production in central Pb + Pb at v^^atat ~ ^T-'i GeV from 
WA9& and in central Au + Au at \/sjvjv = 1^0 GeV, as well as preliminary data at 
200 GeV from PHENI}JSE] and central/peripheral data from STAR E are shown. 
The sum of estimated statistical and systematic errors is indicated. 



February 6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



JET QUENCHING AND RADIATIVE ENERGY LOSSIN DENSE NUCLEAR MATTER31 

Our main results for central Au + Au reactions, which includes three 
important nuclear effects, Cronin+Shadowing+Quenching, are presented in 
Fig. ^1 Jet tomography consists of determining the effective initial gluon 
rapidity density dN^ /dy that best reproduces the quenching pattern of the 
dat a"*^ 5, 24 | gpg ^j^^ large Cronin enhancemenlP^ is reduced by a fac- 
tor of two with dN^ /dy = 350 but the data are more consistent with a 
smaller gluon density < 200. Unfortunately at this low energy the results 
are very sensitive to the details of the model. At RHIC, for > 2 GeV jet 
quenching dominates, but surprisingly the rate of variation with pT of the 
Cronin enhancement and jet quenching conspire to yield an approximately 
constant suppression pattern with magnitude dependent only on the initial 
dN^ /dy. At higher px > 20 GeV the softening of the initial jet spectra due 
to the EMC modification of the PDFs compensates for the reduced energy 
loss. This unexpected interpla y am ong the three nuclear effects at RHIC 
was our main prediction in Reff^. Preliminary -^s = 200 AGeV PHENIX 
and STAR results included in Fig. [THl support the predicted magnitude and 
approximate pr dependence of the suppression. At LHC energies the much 
larger gluon densities dN^ /dy ~ 2000 — 3500 are expected to lead to a 
dramatic variation of quenching with px, as shown. In nuclear media of 
high opacity the mean fractional energy loss (AE) /E of moderately hard 
partons can become of the order of unity. For LHC this may be reflected 
in the px < 10 GeV region through deviations from the extrapolated high- 
er suppression trend. Hadronic fragments coming from energetic jets would 
tend to compensate the rapidly increasing quenching with decreasing trans- 
verse momentum (seen in Fig. I18|l and may restore the hydrodynamic-like 
participant scaling in the soft regime. 



2.6. Enhanced Baryon/Meson Ratios 

One of the unexpected results reported during the first year RHIC run at 
\/^NN ~ GeV was that, in contrast to the strong 7r° quenching for 
2 GeV < pt < 5 GeV, the corresponding charged hadrons were found to 
be suppressed by only a factor ~ 2 — 2.5. Even more surprisingly, the iden- 
tified particle spectra analysis suggests that Rb[pt) = p/tt", p/7r+ > 1 for 
Pt > 2 GeV. Thus, baryon and antibaryon pr oduction may in fact dominate 
the moderate to high px hadron flavor yieldJ^^^SEHEnf These and other 
data point to a possible novel baryon transport dynamics nucleus-nucleus 
reactions. More recent results corroborate the non-perturbative ba ryon p ro- 
duction hypothesis through equally abundant A and A productioi i I I It 
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has also been observed that the mean transverse momentum {pt) b for 
various baryon and antibaryon species is approximately constant and de- 
viates from the common hydrodynamic flow systematics of soft hadron 
production in A + A collisions. Identified particle analyses from the recent 
V^NN — 200 GeVRHIC run find similar puzzling features of moderate-pr 
baryon spec trgl^l^l^. ____ 

In Kefs.™!! the GLV jet energy losjsisni 

was combined 

with a topological non-perturbative baryon production and transport 
mechanism^^ ^ to gain insight into the anomalous anti+baryon behav- 
ior at RHIC. Phenomenological applications are currently based on Regge 
theory where a Regge trajectory J = a{0) + a'( 0)M^ is specified by its 
intercept a(0) and slope a'(0). It has been argued^E^l that aj(0) ~ 0.5 
and a'j{0) ~ 1/3q^(0). Regge theory gives exponential rapidity correla- 
tions, which in the presence of two sources (at ztYmax) lead to net baryon 
rapidity density in central A + A collisions of the form: 

. + «)(! - 0,(0)) ■t,';"-",:'""'' . (52) 

ay smh(l - a,/(0))yniax 

It is evident from Ea. H52(l that the net baryon distribution integrates to 2 A 
and in going to peripheral reactions scales as Npart- At RHIC energies of 
y/s = 130(200) GeV, corresponding to Ymax = 4.8(5.4), in central reactions 
dN^^^/dy = d{p-p)/dy + d{n-n)/dy + d{A-A)/dy + -- ■ ~ 18(13.5). The 
relative contribution of each baryon species can be evaluated from isospin 
symmetry and strangeness conservation via comparison to midrapidity kaon 
production. 

The high px part of the hadron spectra is computed from Eq.©. 
Hadronic transport in small-to-moderate pt is effectively controlled by the 
slope of the Regge trajectory. This would suggest that the baryon/meson 
mean inverse slopes in a phenomenological pT-exponential (~ e'^^^'^) soft 
particle production model are related as Tg : Tm — %/3 : V^. Soft pion 
production, however, is largely dominated by resonance decays, where the 
cumulative effect from the random walk in pr due to string breaking is de- 
stroyed. This leads to the relation (r^rVj. : (pt)k ■ {pt)b — 1 : (l-^\/2) : \/3 
(220 MeV : 275 MeV : 380 MeVpHlHll. One also notes that in the limit of 
pair production dominated by junction-antijunction loops (which we con- 
sider here) the transverse momentum distribution of antibaryons closely 
resembles that of baryons (with {pt)b — {pt)b)- 

We find an enhanced baryon/meson ratio Rb{pt) > 1 that decreases 
with centrality in a finite pt window as illustrated in Fig. 1191 For pT > 
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Fig. 19. The p/n~ ratio versus pj- for 3 centrality classes at \/3jvjv = 1^0 GeV from 
RefsIM A'pari (solid) versus 7Vi,i„ (dashed) geometry is shown. Boosted thermal source 
computation and ratio of fits to PHENIX data are shown for comparison. 



5 GeV the ratio is reduced below unity and approaches the perturbative 
calculation. The centrality dependence of Rb{j>t) is largely driven by the 
non-Abelian energy loss that suppresses perturbative pion production for 
Pt > 2 GeV. It is important to note that an enhanced p/tt ratio is also 
observed in peripheral reactions and is consistent with measurements in 
p-\-p. In Fig. l20l barvon enhancement is reflected in the different suppression 
factor Raa{pt) for neutral pions and inclusive charged hadrons. At high 
transverse momenta the suppression ratio is expected to be similar and 
driven by fragmentation of quenched jets. The rate at which this common 
Raa is approached depends on the fragmentation functions into baryons 
that are currently poorly constrained or unknown. 



2.7. High-pT Azimuthal Asymmetry 

A new way to probe the en ergy loss AE in variable geometries was re- 
cently proposed in Ref.l^^^^. The idea was to exploit the spatial azimuthal 
asymmetry of non-central nuclear collisions. The dependence of AE on the 
path length naturally results in a pattern of azimuthal asymmetry of 
high Pt hadrons which can be measured via the differential elliptic flow 
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parameter (second Fourier coefficient), V2{pt)'. 



V2{pt) 



(pI) + (pD 



(cos2(/)) = 



J^^ d(j) cos 2(j) 



dy pT dpT d(p 



dN*^ 



(53) 



Jo dy pt dpT d(f} 

In Ref.l2£l we predicted V2{pt) for two models of initial condition: 
(specified by dN^ /dy) wliich differ by almost an order of magnitude. The 
longitudinal expansion of the plasma and the corresponding mean energy 
loss were given by Eas. (|39|) . H44|l with a = 1. A novel element of the anal- 
ysis was the discussion of the interplay between the azimuthally asymmet- 
ric soft (hydrodynamic) and hard (quenched jet) components of the final 
hadron distributions. In non-central A + B reactions the low pt hadrons are 
also expected to exhibit azimuthal asymmetry caused by the hydrodynamic 
flow^-'^^'^. We therefore modeled the soft component with the following 
Ansatz: 

dNsib) dns e-PT/^° 



27tT§ 



■ {1 + 2v2s{pt) cos{2(j)) + ■ ■ ■) 



(54) 



dyd'^PT dy 

Here we took Tq « 0.25 GeV and incorporate the hydrodynamic elliptic 
flow predicted in Ref.l^ and found to grow monotonically with pt as 



V2s{pt) ~ tanh(pT/(10 ± 2 GeV)) 



(55) 
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The hydrodynamic elliptic flow was found^to be less sensitive to the initial 
conditions than the high pT jet quenching studied in RefP^. 

We showed that the combined pattern of jet quenching in the single in- 
clusive spectr a a nd the differential elliptic flow at high px provide comple- 
mentary toolJ2^ that can determine the effective density of gluons created 
in the early stages of relativistic heavy ion reactions. 

0.30 
0.25 
0.20 
2 0.15 

CM 

^ 0.10 
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0.00 

1 2 3 4 5 6 
Pt [GeV] 

Fig. 21. The interpolation of «2(pt) between the softhydrodynamiJ^3nH and hard 
pQCD regimes is shown for b = 7 fm adapted from Refl^Sl Solid (dashed) curves corre- 
spond to sharp cylindrical (diffuse Woods-Saxon) geometries. 
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Fig. 12] shows the predicted pattern of high-pT anisotropy. Note the 
difference between sharp cylinder and diffuse Woods-Saxon geometries at 
6 = 7 fm approximating roughly 20-30% central events. While the central 
(6 = 0) inclusive quenching is insensitive to the density profile, non-central 
events clearly exhibit large sensitivity to the actual distribution. In partic- 
ular, the sharp elliptic geometry also simulates one of the possible scenarios 
of large energy loss and predominant surface emissio We conclude that 
V2{pt > 2 GeV, 6) provides essential complementary information about 
the geometry and impact parameter dependence of the initial conditions 
in A + A reactions and the magnitude of jet quenching. In particular, the 
rate at which the f 2 coefficient decreases at high px is an indicator of the 
diffuseness of that geometry. Minimum bias STAR data at RHIC^ now 
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confirm the predicted saturation of V2{pt)- For pr > Q GeV data is still 
inconclusive due to the large error bars and can still accommodate different 
scenarios: constant versus slightly decreasing V2{pt)- High statistics mea- 
surements in the future RHIC runs can help t o experimentally resolve this 
question and further test the predicted""'"^ slow decrease of V2 at large 
transverse momen ta. 

hi JRe fs N^ l ^^i in an independent analysis the GLV radiative energy 
losJ^^EHl of fast partons was coupled to a dynamically evolving soft hy- 
drodynamic background to obtain a quantitative description of the disap- 
pearance of the back-to-back jet correlation^ 

3. The WW Approach: Parton Energy Loss with Detailed 
Balance 

The above studies of parton energy loss have concentrated on gluon radia- 
tion induced by multiple scattering in a medium. Since gluons are bosons, 
there should also be stimulated gluon emission and absorption by the propa- 
gating parton because of the presence of thermal gluons in the hot medium. 
Such detailed balance is crucial for parton thermalization and should also 
be imp orta nt for calculating the energy loss of an energetic parton in a hot 
mediuirP^. 

3.1. Final-state Absorption 

Let us assume that the hot medium is in thermal equilibrium shortly after 
the production of the hard parton. One should then take into account of 
both stimulated emission and thermal absorption in the final state radiation 
of a jet in a thermal medium with finite temperature T. One has then the 
probability of gluon radiation with energy w, 

Ng{zE)S{Lo + zE) 

+ (1 + Ng{zE)) 6{oj - zE)0(l - z)] P(|) , (56) 

where, 7Vg(|k|) = l/[exp(|k|/T) — 1] is the thermal gluon distribution. We 
have define the splitting function Pgq{z) = P{z)/z = [1 -|- (1 — z)'^]/z for 
q — !■ gq. The first term is from thermal absorption and the second term from 
gluon emission with the Bose-Einstein enhancement factor. For E ^ T, one 
can neglect the quantum statistical effect for the leading parton. Note that 
the vacuum part has a logarithmic infrared divergence while the finite- 
temperature part has a linear divergence, since 7Vg(|k|) ~T/|k| as |k| ^0. 
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These infrared divergences will be canceled by the virtual corrections which 
also contain a zero-temperature and a finite-temperature part. In addition, 
the virtual corrections are also essential to ensure unitarity and momentum 

conservation in the QCD evolution of the fragmentation functions. However, 
they do not contribute to the effective parton energy loss. The remaining 
collinear divergence in the above spectrum can be absorbed into a renor- 
malized fragmentation function that follows the QCD evolution equations. 

Subtracting the gluon radiation spectrum in the vacuum, one then ob- 
tains the energy loss due to final-state absorption and stimulated emission, 



du)u) 



asCp 
27r 



T=0 



E / dz 



ki 



+P{z)Ng{zE)6{l - z) 



P{-z)Ng{zE) 



(57) 



Even though the stimulated emission cancels part of the contribution from 
absorption, the net medium effect without rcscattering is still dominated 
by the final-state thermal absorption, resulting in a net energy gain, i.e. a 
negative energy loss. For asymptotically large parton energy, E ^ T, one 
can complete the above integration approximately and have, 



abs 

E 



3 E^ 



In- 



iET 



■ 2 - 7e + 



6C'(2) 



(58) 



where, 7e ~ 0.5772 and C'(2) ~ —0.9376. The quadratic temperature de- 
pendence of the leading contribution is a direct consequence of the par- 
tial cancellation between stimulated emission and thermal absorption, each 
having a leading contribution linear in T. 



3.2. Induced Absorption 

As in the case of final-state absorption, one can also include stimulated 
emission and thermal absorption when calculating the induced radiation 
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probability at the first order in opacity, 



duj 



Co 



SndAdii. 
N 



dz 

z 



Tr 



(27r)2 
f 2i?e 



(27r)2 



(1 + Ng{zE)) 5{lu - zE)0{l ~ z) + Ng{zE)S{uj + zE) 



asC2CFCA 



dz 

z 



dkl 



d'^q± 
(2n 



N 



(i?e(l - e^'^i^")) Ng{zE)5{uj + zE) 



+ (1 + Ng{zE)) 5{uj - zE)e{l - z) 



(59) 



The factor 1 — exp(iCi;ij/io) reflects the destructive interference arising from 
the non-Abelian LPM effect. The target distribution is assumed to be an 
exponential form p{y) = 2 exp{—2y / L) / L. 

As in the final-state absorption, the contribution from thermal absorp- 
tion associated with rescattering is larger than that of stimulated emission, 
resulting in a net energy gain. However, the zero-temperature contribution 
corresponds to the radiation induced by rescattering which will lead to an 
effective energy loss by the leading parton. This is the energy loss obtained 
the previous sections and we denote this part as AE^^J^. The remainder or 
temperature-dependent part of energy loss induced by rescattering at the 
first order in opacity is then defined as 



ae[11 = I doj 



dp^^) dP^^^ 



dhJ 



djjj 



T=0 



(60) 



which mainly comes from thermal absorption with partial cancellation by 
stimulated emission in the medium. According to Ea. H59|) . 
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rad 
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(61) 
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- P(-z)(i?e(l-e"^i'^i°))], 



(62) 
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In the limit q±max ~^ oo, the angular integral can be carried out by partial 
integration. These contributions to the energy loss become 
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where, 7 = ii^L/{AzE) and 
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^(7) 
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dzNg{zE)h{-f) \p{z)e{l - z)- P{~z) 

:[§ - arcsin(27)] , 7 < 1/2 



i2= ln[27 + v/4^23T] , 7 > 1/2 . 



(63) 
(64) 

(65) 



One can approximate h{j) with tt^ + (11/4 — 27r)7^ + (5/2)7"^ for 7 < 1/2 
and ln(47) + O.I/7 + 0.028/7^ for 7 > 1/2. In the limit of £;L > 1 and 
3> ^, One can then get the approximate asymptotic behavior of the 
energy loss, 



AE, 



,(1) 

rad 
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2_E/ 

In ^r- - 0.048 



(1) 
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3 XgE^ 



In 



7e 



6C'(2) 



(66) 
(67) 



Again, the thermal absorption results in an energy gain ( or negative energy 
loss). This result is accurate through the order of 1/_E. In Ea. H67|) . we have 
assumed is^L/T ^ 1 and kept only the first two leading terms. In this limit, 
the average formation time for stimulated emission or thermal absorption is 
much smaller than the total propagation length. Therefore, the energy gain, 
Arabs' by thermal absorption (with partial cancellation by the stimulated 
emission) is linear in L, as compared to the quadratic dependence in the 
zero-temperature case. However, the logarithmic dependence on fi^L/T, as 
compared to the factor h\{AET / ji^) in Ea. H58|l for no rescattering, is still a 
consequence of the LPM interference in medium. A quadratic i-dependence 
of Ai;^2 will arise when /x^L/T < 1. 

To evaluate the effect of thermal absorption numerically, we assume the 
Debye screening mass to be /^^ = AnagT"^ from the perturbative QCD at 
finite temperaturJSSI The mean-free-path for a gluon A„ in the GW model 
is72. 



(68) 
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where C(3) ~ 1.202. With fixed values of L/\g and a^, AE'Z/i should be a 
function oi E/n only. Shown in Fig. [22 are ratios of the calculated radiative 
energy loss with and without stimulated emission and thermal absorption 
as functions of E/ fi for L/Xg = 3,5 and — 0.3. Shown in the inserted box 
are the energy gain via gluon absorption with (AE'^J^^) and without {AE^^J^) 
rescattering. For partons with very high energy the effect of the gluon ab- 
sorption is small and can be neglected. However, the thermal absorption 
reduces the effective parton energy loss by about 30-10% for intermediate 
values of parton energy. This will increase the energy dependence of the ef- 
fective parton energy loss in the intermediate energy region. The observed 
flat or slightly px dependence of the hadron spectra suppression indicates 
a strong energy dependence of the parton energy loss. This might indicate 
the effect of the thermal absorption, in particular in the intermediate pt 
region. 
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4. The WOGZ Approach: Parton Energy Loss and 
Modified Fragmentation Functions in Nuclei 

So far we have reviewed the parton energy loss of an energetic parton in- 
duced by multiple scattering in a hot QCD medium. The effective parton 
energy loss is shown to be proportional to the gluon density. Therefore mea- 
surements of the parton energy loss will enable one to extract the initial 
gluon density of the produced hot medium within our assumed screened 
potential model for multiple parton scattering in a dense gluonic medium. 
Such an exercise will be more robust if one can also measure the parton 
energy loss and the gluon density inside a cold nuclear medium. For this 
purpose, one has to rely on other complementary experimental measure- 
ments such as parton energy loss in deeply inelastic scattering (DIS) off 
nuclear targets. One can then study parton energy loss and extract the 
initial gluon dens ity in high-energy heavy-ion collisions relative to that in 
a cold nucleuJ^. 

Extending the screened potential model of multiple parton scattering to 
a cold nuclear medium is somewhat problematic because it applies mainly 
to a medium where colors are screened but not confined. In this section 
we will study parton multiple scattering inside a nucleus where colors are 
confined to the size of a nucleon within generalized factorization theory 
of pQCD. In this case, parton propagation will be different from that in 
a partonic medium and one should expect the parton energy loss to be 
related to the nucleon size or confinement scale which, as we shall show, 
can be related to the gluon density inside a cold nucleus. In addition, one 
can calculate directly the modification of the fragmentation functions due to 
multiple scattering and induced gluonbremsstrahlung. This section reviews 
developments published in RefsllSSlS^E^^S^ which are referred to here as 
WOGZ and are applied in RefP to cold and hot nuclei. 

In contrast to the situation in QED, the energy loss of a parton cannot 
be directly measured because partons are not the final experimentally ob- 
served particles. The total energy of a jet, traditionally defined as a cluster 
of hadrons in the phase space, will not change much due to medium in- 
duced radiation, because a jet so defined contains particles both from the 
leading parton and from the radiated gluons. This is particularly the case 
if multiple scattering and induced radiation do not dramatically change 
the energy profile of the jet in phase-space. It is also virtually impossible 
to determine the jet energy event by event because of the large background 
and its fluctuation in heavy-ion collisions. One then has to resort to parti- 
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cle distributions within a jet and study the effect of parton energy loss by 
measuring the modification of the fragmentation function of the produced 
parton, Da^h{z, jJ?), where z is the fractional energy of the parton a carried 
by the produced particles h. 

Since the produced quark in DIS is far off-shell, the final-state radia- 
tion leads to the scale dependence of the fragmentation functions as given 
by Dokshitzcr-Gribov-Lipatov-Altarelli-Parisi (DGLApjinni QCD evolu- 
tion equations. When a parton is produced in a medium, it will suffer 
multiple scattering and induced radiation that will give rise to an addi- 
tional term in the DGLAP evolution equations. This then leads to the 
medium modification of the DGLAP evolution of the parton fragmentation 
functions. As a consequence, the modified fragmentation functions become 
softer. This can be directly translated into the energy loss of the leading 
quark. 



4.1. Generalized Factorization 

To study parton energy loss in eA DIS, we consider the semi-inclusive pro- 
cesses, e(Li) -I- A{p) — > 6(^2) + h{ih) + X, where Li and L2 are the four- 
momenta of the incoming and the outgoing leptons, and £h is the observed 
hadron momentum. The differential cross section for the semi-inclusive pro- 
cess can be expressed as 

^^^^'"d^Z^ - ^Qi^M.^^.^ . (69) 

where p ^ [p'^,0,Oj_] is the momentum per nucleon in the nucleus, q ~ 
L2 — £1 = [—Q'^/2q^,q^,0j_] the momentum transfer, s ^ (p + LiY and 
ttEM is the electromagnetic (EM) coupling constant. The Icptonic tensor is 
given by L^^/ — l/2Tr(7 • ^17^7 • L2^u) while the semi-inclusive hadronic 
tensor is defined as, 

i?^.^ = lY.(^\M0)\x,h){x,h\M0)\A) 

X 2TTS\q + p ~ px - ih) (70) 

where J2x runs over all possible final states and — J2q ^q'4'ql(i.'4'q is the 
hadronic EM current. 

In the parton model with the coUinear factorization approximation, the 
leading-twist contribution to the semi-inclusive cross section can be factor- 
ized into a product of parton distributions, parton fragmentation functions 
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and the partonic cross section. Including all leading log radiative correc- 
tions, the lowest order contribution (0(0")) from a single hard 7* +(? scat- 
tering can be written as 



dzh 

9 



Y^el J dxf^ix,^^j)Hj^J{x,p,q)Dg^hizh,^^i^); (71) 



1 27r 

where the momentum fraction carried by the hadron is defined as Zh = 
and xb = Q^/2p'^q~ is the Bjorken variable, fij and are the 
factorization scales for the initial quark distributions f^{x,iij) in a nu- 
cleus and the fragmentation functions Dq^hizh, fJ-^), respectively. Consid- 
ering the leading logarithm approximation of the radiative correction to the 
fragmentation process as shown in Fig. 1231 the renor mali zed quark frag- 
mentation function Z3g^/i(z^, /i^) satisfies the DGLAF^^ QCD evolution 
equations, 

^ln^^^ = J^^ - K^,,(.)i?,^.(../.,M ) 

+ 1q^gq{z)Dg^h{zh/z,fi^)] , (73) 



where 



F 



1 + 2^ 3 



(74) 



(l-z)+ 2 

lq^gq{z) = 79^95 (1 " z) (75) 

are the splitting functions. The '+' function is defined as 



(1 - ^)+ 

with F{z) being any function which is sufficiently smooth at z = 1. 

In a nuclear medium, the propagating quark in DIS will experience 
additional scatterings with other partons from the nucleus. The rescatter- 
ings may induce additional gluon radiation and cause the leading quark to 
lose energy. Such induced gluon radiations will effectively give rise to ad- 
ditional terms in the evolution equation leading to the modification of the 
fragmentation functions in a medium. These are the so-called higher-twist 
corrections since they involve higher-twist parton matrix elements and are 
power-suppressed. We will consider those contributions that involve two- 
parton correlations from two different nucleons inside the nucleus. They 
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xp 
Ap 



xp 



Ap 



Fig. 23. The hard partonic part of the first order radiative correction to the fragmen- 
tation process. 



are proportional to the size of the nucleuff^ and thus are enhanced by a 
nuclear factor A^^^ as com pared to two-parton correlations in a nucleon. As 
m previous 

studiej™, 

we will neglect those contributions that are not 
enhanced by the nuclear medium. 




Fig. 24. A typical diagram for quark-gluon re-scattering processes with three possible 
cuts, central(C), left{L) and right(R). 



We will employ the generalized factorization of multiple scattering 
processes . In this approximation, the double scattering contribution to 
radiative correction from processes like the one illustrated in Fig. can 
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be written in the following form, 



dW, 



D 



dzh 



E 

9 



1 dz 



-D 



q — >h 



dy 

2tt 



dVi dy. 



2 (2^)2 



d^kj 



X e 



-ikr-yT 



-rrD 



H^uiy ,yi ,y2 ,kT,p,q,z); 



X -{A\%{0)j+A+iy-,OT)A+iy-,yT)My"M) ■ (77) 



Here H^^{y ^y^ , j/j : ^TtP, Q, z) is the Fourier transform of the partonic 
hard part H^i,{x, Xi,X2, kx^p, q, z) in momentum space, 



Hn,,iy ,yi ,y2 ,kT,P,q,z) = / dx 



dxi dx2 



^ixip^y +ix2P y-i 

1-K 27r " 

xe*(^-^^--^)f"^^"i?^,(x,xi,a;2,A:T,p,g,z) , (78) 

where kx is the relative transverse momentum carried by the second parton 
in the double scattering. Values of the momentum fractions a;, x\ and x^ 
are fixed by 5- functions and poles in the partonic hard part. They normally 
depend on fc^. 

In order to pick up the next-leading-twist contribution, we expand the 
partonic hard part around fcy = 0, 

H^Ay^ ,yi ,y2 ^f^T:P,q,z) = H^^y- ,y^ ,y:^ ,kT = 0,p,g,z) 



dH 



■D 



dkS 



1 dm: 



kT=0 



2 9fcf a4 



(79) 



This is known as the coUinear expansior^^. On the right-hand-side of 
Ea. (|79|l . the first term gives the eikonal contribution to the leading- twist 
results. It does not correspond to the physical double scattering, but simply 
makes the matrix element in a single scattering gauge invariant. The second 
term for unpolarized initial and final states vanishes after being integrated 
over kx- The third term will give a finite contribution to the double scatter- 
ing process. Substituting Ea. (|7^ into Ea. ((77|) . integrating over d^kx and 
d^yx, we obtain 



dW, 



dzh 



E 

1 



1 dz 



— Dq^h{zh/z) / -^dy^ dy. 



dy 



X - (A| (0) 7+ F+ {y^ ) F+'^ (y^ ) (y" ) | A) 



2 dk!f,dk^ 



,(80) 



kT=0 
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where k^A+k!^A+ are converted into field strength F'^+p0+ by partial 
integrations. 

4.2. Double Parton Scattering 

There are all together 23 cut-diagrams that contribute to the leading twist- 
four corrections to the quark fragmentation function in eA DIS. For simpli- 
fication of the calculation, one can use the helicity amplitude approximation 
by Guo and Wang-SS^'^in the limit of soft gluon radiation. Such an approx- 
imation enables one to simplify the calculation of the radiation amplitudes. 
However, a complete calc ulation of the all cut-diagrams were recently done 
by Zhang and WanJ^^^. Before we list the results, we first consider the 
contribution from Fig. El in detail. 

Using the conventional Feynman rule, one c an wri te down the hard 
partonic part of the central cut-diagram of Fig. |2 jP'^ l ^^^ l 

Hc,.iy-,yi,y2,kT,P,q,z)= f dx^^e^^^'^^y'+^^^p^y^ 

J ZTT ZTT 

xe'^^-^'-^'^P^y^ f S^Itt \p ■ ^-/^p'^pPH^p^, 



(27r)4 2 

x2Tr6+{e^) S{1 - z - —) . (81) 

q 

fr _ Cf_ 4 7 ■ (g + Xip) -1 ■ {q + xip- £) 

2Nc (q + xip)'^ — le (g + xip — — le 

^= W(, + ,p_^)2+,,7/.(, + ,p)2+,,2.^+(^,)' (82) 

where e"'^(£) is the polarization tensor of a gluon propagator in an axial 
gauge, n-A = with n ~ [1, 0^, 0^], and £q — q+{xi+X2)p+kT—i are the 
4-momenta carried by the gluon and the final quark, respectively, z = (~ /q~ 
is the fraction of longitudinal momentum (the large minus component) 
carried by the final quark. 

To simplify the calculation, we apply the coUinear approximation to 
complete the trace of the product of 7-matrices, 

p-H^ppP^j.iq-LTJj-p'^H^ppp] . (83) 

After carrying out momentum integrations in x, xi, X2 and with the 
help of contour integration and (5-functions, the partonic hard part can 
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be factorized into the product of 7-quark scattering matrix Hjj^J {x,p, q) 
[Eq.O] and the quark-gluon rescattering part , 

-ttD 



Hf,^{y ,y2 ,kT,p,q,z) ^ J dxHj^J {x,p,q) 

^ H° {y^ ,Vi ,V2 ikT,x,p,q,z) . (84) 

Contributions from all the diagrams have this factorized from. Therefore, 
we will only list the rescattering part H for differen t diag rams in the 
following. For the central-cut diagram in Fig. [21 it readi 



Hc{Fig\I^{.v ,yi ,y2 ,kT,x,p,q,z) ^ j -jr^^F ^ _ ^ 

2iias— _ _ _ 

x^r^-^CffM[^(2/ ,yi ,y2 ^^T,kT,x,p,q,z) , (85) 

^d{y' - yr)(i - e-"'-p*yi){i - e-'''^p'^^y'^yi'>) . (86) 

Here, the fractional momentum is defined as 

4 fcr ^ 2kT ■ It 

2p+q z[l — z) 1 ^ 

and x = Xb = I'^V^I^ is the Bjorken variable. 

The above contribution resembles the cross section of a dipole scattering 
and contains essentially four terms. The first diagonal term corresponds to 
the so-called hard-soft process where the gluon radiation is induced by 
the hard scattering between the virtual photon and an initial quark with 
momentum fraction x. The quark is knocked off-shell by the virtual photon 
and becomes on-shell again after radiating a gluon. Afterwards the on-shell 
quark (or the radiated gluon) will have a secondary scattering with another 
soft gluon from the nucleus. The second diagonal term is due to the double 
hard process where the quark is on-shell after the first hard scattering with 
the virtual photon. The gluon radiation is then induced by the scattering of 
the quark with another gluon that carries finite momentum fraction Xh+xo. 
The other two off-diagonal terms are interferences between hard-soft and 
double hard processes. In the limit of collinear radiation (xj^ — > 0) or when 
the formation time of the gluon radiation, r/ = Xjx^P^ is much larger than 
the nuclear size, the two processes have destructive interference, leading to 
the LPM interference effect. 
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One can similarly obtain the rescattering part H of other central-cut 



diagrams (a-d) in Fig. 123 



Hc{,a){y ,yi ,y2 ^>^T,x,p,q,z) = J ■ 



{£t - kT)2 27r 1-z 

X ^^^C(a) {y^ ,yi,y2i^T,kT,X, p, q, z) , 



Ida) (y- , yr ,y2,^T,kT,x,p, q, z) = e^i-+--)p^y- (y^ -y^^9{~y-) 



Hci,) {y-,y^,y2,kr,x,p,q,z)^J _ ^ J^^^^^^, 2^ 



27ras— , _ _ _ 
'x^^Ic{b)[V ,yi ,y2 JT,kT,x,p,q,z) , 

7cmiy~.yi,y2.f^T,kT,x,p,q,z)^e^(-+--^p^y'+^--p^(y^-y^^e{~y^) 

>^0{y- - y-)e--'^^P^^y'-yi') "^P^y^ , (89) 
jjD _ _ f {It ~ kT) ■ (^T - (1 ^ ^)kT) 

Hc^y ,y,.y2. kr. x,p^ q,z) = J _ 

Us Ca 1 + z'^ 2Tras-r ^ 
~ 1 _ ^ -J^^c{c)[y ,Vi,V2 J^T,kT,x,p,q,z), 

xe{y- -y-)e-''=^P*y2 



T7^ . - - - J. X f ,,2 (^T - kT) ■ (^T - (1 - ^)kT) 

i^c(.)(y ,yi,y2. kr, x,p, q,z)^ I _ ^^^^^^^ ^^^^.^^ 



Us Ca1 + z^ 2TTas- . 

ic(d)[y ,yi ,y2 ,>^T,kT,x,p,q,z), 



2tt 2 1-z N, 



ic(d){y-,yi,y2,eT,kT,x,p,q,z) - e'(-+--)f y'+^-^p ivi-y-,) e(^-y-) 
xe{y- ~y^)e-''=^p*^y'-y^'> 



To complete the calculation we also have to consider the asymmetrical- 
cut diagrams (left-cut and right-cut) that represent interferences between 
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single and triple scatterings. They can be obtained with similar procedures. 
We refer readers to Ref."'^'-'^ for a 1 
those asymmetrical-cut diagrams. 



We refer readers to Ref."'^'-'^ for a list of the the rescattering part of all 
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Fig. 25. Four central-cut diagrams that contribute to the final results. 



To obtain the double scattering contribution to the semi-inclusive pro- 
cesses of hadron production in Eq. H8Q(I , one will then have to calculate the 
second derivatives of the rescattering part H . After a closer examination 
of these rescattering parts, one can find that all contributions from the 
asymmetrical-cut diagrams have the form 

(t ■ {It - /(^)kT) 



H 



■D 

asym 



4( 



iXp+Y- 



(92) 



where f{z) — 0, 1, 1 — z, z, X is the longitudinal momentum fraction 
and represents the spatial coordinates. One can prove that the second 
derivative of the above expression vanishes at fc^ = 0, 

„2 ^T-(^T-/WkT) 



4(^T - /(2)kT)2 



= 0. 



(93) 
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Therefore, all contributions from the asymmetrical-cut (right-cut and left- 
cut) diagrams will vanish after we take the second partial derivative with 
respect to kx when we keep only the leading terms up to 0{xb /Q'^^t)' 



(94) 



In the same way we find that some of the central-cut diagrams will not 
contribute to the final results, either. In fact, after taking the second partial 
derivative with respect to kr only four central-cut diagrams shown in Fig. 1251 
will contribute to the final result. 

Including only those contributions that do not vanish after the second 
derivative with respect to kx, we have 



27r 1 - z 
'^'^^(l_e-"^P+!^2")(l-e 



^(-2/2")%"-yr) 

) 



4Cf(1 z) ixLp+(y -y^ )p-iXLP'^y2 



+ 



2Ca(1-z) 



g-ixip+i/j _ ^-ixLP^{y -y^ ) 



2Ca(1 ^') ^-ixLp^(y -y^ _ ^-ixLP+y^ ) 



+ 0{xB/Q^i^T)] 



(95) 



The first term at the right-hand side in Ea. H95|l comes from the contri- 
bution of ^c(a) which is the main contribution in the helicity ampli- 
tude approximatioJ^^I^. It contains hard-soft, double hard processes and 
their interferences. The other three terms come from diagram (b),(c),(d) of 
Fig. 123 respectively. They constitute corrections to the first term in powers 
of 1 — z. The second term that is proportional to (1 — z)'^ is from the final 
state radiation from the quark in the double hard process in Fig. l25f bl The 
third and fourth terms are the results of the interference of the final state 
radiation from the quark and other radiation processes (initial state radi- 
ation and radiation from the gluon line). They contain both double hard 
processes and interferences between hard-soft and double hard processes in 
Fig.EHc) and (d). 

Substituting Ea. H95|l into Eas. 1)84(1 and ((8()|l . we have the semi-inclusive 
tensor from double quark-gluon scattering including the contribution be- 
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yond the helicity amplitude approximation, 

"1 dz 



dzh 



Y: Jd.H(^Ji.p,q)J^ ^D,^,iz,/z)^CA\ 



1 + z' 



di'^ 27ras 



T,^,(x,XL) + (l-z)T^fi)(a;,a;L) 



+ ^il-z)XTi-,x,) 



where 



Tqg{x,XL) 



dy 
27r 



(96) 



X e^^^+^-^P^y-e{-y^)d{y- -y^) 



(97) 



rpA{l 



dy' 



qg \X,Xl) = I -l^'^Vl '^y2 ^ 



i{x+XL)p y 



2g-^2;ip+(j/ -y-^ +y^ ) 



1 



{-V2 Wiv - Vi ) 



X - {A\i^q (0) 7+ F+ {y^ ) (yf )^q{y-)\A) 



(98) 



T^}'^\x,Xl)= I — dy-dy-e'^'P^y +^^lp+ {y^ -y^ )i 

^y J Itx 



~V2 )^(y - vi ) 



(99) 



are twist-four parton matrix elements of the nucleus. Evidently these par- 
ton matrix elements are not independent of each other. T^{x, xl) has the 
complete four terms of soft-hard, double hard processes and their inter- 
ferences. Therefore it contains essentially four independent parton matrix 
elements. r^*-^-'(x, x^) and T^^'^\x, xl) are the results of the corrections 
beyond the helicity amplitude approximation. But these two matrix ele- 
ments are already contained in r^(a;, xl)- 

During the collinear expansion, we have kept £t finite and took the limit 
kr ^ 0. As a consequence, the gluon field in one of the twist-four parton 
matrix elements in Eas. H97|l - 1)99(1 carries zero momentum in the soft-hard 
process. However, the gl uon di stribution xfg{x) at a; = is not defined in 
QCD. As argued in Refs.l^^^^, this is due to the omission of higher order 
terms in the collinear expansion. As a remedy to the problem, a subset 
of the higher-twist terms in the collinear expansion can be resummed to 
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restore the phase factors such as exp(ia;Tp+7y^), where xt = (k^)/2p~^q~z 
is related to the intrinsic transverse momentum of the initial partons. As a 
result, soft gluon fields in the parton matrix elements will carry a fractional 

momentum xt- 

Using the factorization 

approximatioJHl7U|102|101| 

we can relate the 

twist-four parton matrix elements of the nucleus to the twist-two parton 
distributions of nucleons and the nucleus, 

_C 

XA 

+ f^{x)ixL + XT)f^{xL + xt)] , (100) 

where C is a constant, xa = ^/MRa, fq{x) is the quark distribution 
inside a nucleus, and f^{x) is the gluon distribution inside a nucleon. A 
Gaussian distribution in the light-cone coordinates was assumed for the 
nuclear distribution, p{y^) — poexp{y~ /2i?^ ), where i?^ = \j2RaMIp^ 
and M is the nucleon mass. We should emphasize that the t wist-fo ur matrix 
element is proportional to XjxA = RaM, or the nuclear siz m. 

Notice that the off-diagonal matrix elements that correspond to the 
interferences between hard-soft and double hard processes are suppressed 
by a factor of exp{—x'j^/x\). This is because in the interferences between 
double-hard and hard-soft processes, there is actually momentum flow of 
xlP~^ between the two nucleons that the initial quark and gluon come from. 
Without strong long range two-nucleon correlation inside a nucleus, the 
amount of momentum flow x^p^^ should then be restricted to the amount 
allowed by the uncertainty principle, 1 / ^ p^ / RaM . Similarly, the other 
two-parton matrix elements in Eas. (|98() and (|99|l can be approximated as 

+ f^{x)ixL + XT)f^{xL + xt)] e-^^/^- 

- 2f^{x){xL + XT)f^{xL +xt)} , (101) 

T^f'Hx,XL) - —f^(x){xL + XT)f^{xL + xt) ■ (102) 

Xa 

^From the above estimate of the matrix elements, both T^g{x, xl) and 
T^^^\x, xl) contain a factor 1 — e~^i^^^^ because of the LPM interference 
effect. Such an interference factor will effectively cut off the integration over 
the transverse momentum at xl ^ xa in Ea. H96|l . As we will show later 
in the calculation of the effective energy loss, the integration with such a 
restriction in the transverse momentum due to LPM interference effect will 



February 6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



JET QUENCHING AND RADIATIVE ENERGY LOSSIN DENSE NUCLEAR MATTER53 



give rise to a factor 1/xa in addition to the coefficient f^{x)/xA- Conse- 
quently, contributions from double scattering in Eq. that are associated 
with T^g{x,XL) and Tqg^\x,XL) will be proportional to R\f^{x). These 
are the leading double scattering contributions in the limit of a large nu- 
clear size. On the other hand, the third term Tqg'^\x^ xl) in Ea. H96|) . which 
does not contain any interference effect, will only contribute to a correction 
that is proportional to RAf^{x). In the limit of a large nucleus, A^^^ ^ 1, 
we will neglect this term in our study of the double scattering processes. 



4.3. Virtual Corrections 

So far we have not considered virtual corrections which will ensure the final 
result to be infrared safe. The calculation of the virtual corrections can be 
calculated similarly as the radiative corrections. On the other hand, they 
can also be obtained via unitarity requirement. When cast into the DGLAP 
evolution equation as Eq. H73() , the real corrections can be interpreted as the 
probability for the quark to radiate a gluon with momentum fraction 1 — z. 
Then one must also take into account the probability of no gluon radiation 
in the evolution to ensure unitarity. Such unitarity requirement gives rise 
to the same virtual correction as calculated from the virtual diagrams. 
The virtual contribution to the quark fragmentation in double scattering 
processes is, for example. 



dzh 



J2 J dxHjS{xp,q)Dq^^ (Zh) 



where 



T^f EE Tf^{x,XL) + (1 ~ z)T^^^<^'\x,xl) . (104) 



Assuming a function F{z) which is sufficiently smooth at z = 1, One can 
single out the infrared divergent part of the following integral, 

i + z^ 2 

/ dz^—F{z)^F{l) dz Ai^; 

1 

/ dz- [2F(l)-(l + z2)F(z)] . (105) 

Jo 1 — 2 

The second term A_F is finite since F{z) is a smooth function of z. 
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We can apply this procedure to the integral in the modified fragmen- 
tation function. The divergent term can be combined with the radiative 
contribution in Eq.JHEl to cancel the infrared divergency. With the help of 
the '+'functiortl^, the final result can be expressed as 



dzh 



E 

9 



dxHl^J{xp,q) 



dl^ dz I 



1 + 



19 



(106) 



dz 





-{l + z^)T^}-^\x,xA . (107) 



Here the implicit z-dependence of T^J'™\x, xl) plays an important role in 
the final result. The above integrand will be proportional to the splitting 
function for a single scattering if one ignores the z dependence of T^(a;, xl )- 
Similarly, the final result for contributions from gluon fragmentation is 

dzh 



^ + " '")(x, XL) + 5{z)^T^}^\x, l\) 



2+ 



, (108) 



where we have used the fact that xl in Ea. H87|) is invariant under the 
transform z —* \ — z and so is T^JyX^ xl)- 



4.4. Modified Fragmentation Function and Parton Energy 
Loss 

Including these virtual corrections and the single scattering contribution, we 
can rewrite the semi-inclusive tensor in terms of a modified fragmentation 
function Dg^h{zh, ^J.'^), 



'1- 
dzh 



^Y^j dxff{x,fif)Hj^J{x,p,q)D,^h{zh,fi') (109) 



where f^{x, /ij) is the quark distribution function which in principle should 
also include the higher- twist contributioi^^ of the initial state scattering. 
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The modified effective quark fragmentation function is defined as 

P d^rp OLq F dZ r . , rj , , . . 

+ / -JT^ —[A-/g^gg{z,X,XLjT)Dq^h[Zh/z) 
Jo *-T Jzh ^ 

+ A'yq^gg{z,X,XL,iT)Dg^h{zh/z)] , (110) 



where Dq^hizh, ^J^) and Dg—,h{zh-:^J^) are the leading-twist fragmentation 
functions. The modified sphtting functions are given as 

1 + 2^ 

l^lq^qg{z,X,XLAT) 



T^}^\x,XL)+5{l-z)AT^}-^\x,ll) 
2'KasCA 



(l-z)+ 



(111) 



Mq^gqiz^x^XLAx) = A7q^qg(l - z,x,xlAt)- (112) 

To further simphfy the calculation, we assume <C <C x. The 
modified parton matrix elements can be approximated by 



T^}"'\x,xl) « —(1 - e~^'^/^'-)f^{x) 

XA 



1 

1 



(113) 



where C = 2CxTf^{xT) is a coefficient which should in principle depend 



' g 

-)2 



on Q and xt- Here we will simply take it as a constant. 

In the above matrix element, one can identify 1/ x^p^ — 2q^ z{\ — z)/ ^"^ 
as the formation time of the emitted gluons. For large formation time as 
compared to the nuclear size, the above matrix element vanishes, demon- 
strating a typical LPM interference effect. This is because the emitted gluon 
(with long formation time) and the leading quark are still a coherent system 
when they propagate through the nucleus. Additional scattering will not 
induce more gluon radiation, thus limiting the energy loss of the leading 
quark. 

The reduction of phase space available for gluon radiation due to the 
LPM interference effect is critical for applying th e LQ S formalism to the 
problem in this paper. In the original LQS approacfP^, the generalized fac- 
torization for processes with a large final transverse momentum l'^ ~ 
allows one to consider the leading contribution in 1/Q^, which is enhanced 
by the nuclear size Ra ~ A^^^. For large and A, the higher- twist con- 
tribution from double parton rescattering that is proportional to asRA/Q"^ 
will then be the leading nuclear correction. One can neglect contributions 
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from more than two parton rescattering. In deriving the modified frag- 
mentation functions, we however have to take the leading logarithmic ap- 
proximation in the limit i'^ ^ Q^, where £t is the transverse momentum 
of the radiated gluon. Since the LPM interference suppresses gluon radi- 
ation whose formation time (r/ ~ Q'^ /^tP^) larger than the nuclear 
size MRa/p~^ in our chosen frame, i'^ should then have a minimum value 
of i"^ - Q^/MRa Q^/A^/^. Here M is the nucleon mass. Therefore, 
the logarithmic approximation is still valid for large nuclei [MRa ^ !)■ 
In the meantime, the leading higher-twist contribution proportional to 
asRA/^T ^ oisR\/Q^ can still be small for large so that one can 
neglect processes with more than two parton rescattering. The parame- 
ter for twist expansion of the fragmentation processes inside a nucleus is 
thus UsA^^^ IQ'^ as compared to UgA'^/'^ /Q'^ for processes with large final 
transverse momentum as studied by LQS . This is why the nuclear mod- 
ification to the fragmentation function as derived in this paper depends 
quadratically on the nuclear size Ra- 

Because of momentum conservation, the fractional momentum in a nu- 
cleon is limited to xl < 1- Though the Fermi motion effect in a nucleus 
can allow o^l > 1, the parton distribution in this region is still significantly 
suppressed. It therefore provides a natural cut-off for in the integration 
over z and £t in Ea. Hll()() . With the assumption of the factorized form of 
the twist-4 nuclear parton matrices, there is only one free parameter C{Q^) 
which represents quark-gluon correlation strength inside nuclei. Once it is 
fixed, one can predict the z, energy and nuclear dependence of the medium 
modification of the fragmentation function. Shown in Figs. IJEl and 123 are 
the calculated nuclear modification factor of the fragmentation functions for 
i^iV and ^^Kr targets as compared to the recent HERMES data^^S, There 
are strong correlations among values of Q^, v and z in the HERMES data 
which are also taken in account in our calculation. The predicted shape of 
the z- and i^-dependence agrees well with the experimental data. A remark- 
able feature of the prediction is the quadratic A'^^^ nuclear size dependence, 
which is verified for the first time by an experiment. This quadratic depen- 
dence comes from the combination of the QCD radiation spectrum and the 
modification of the available phase space in ix or xl due to the LPM inter- 
ferences. Note that the numerical results shown here are obtained with the 
original helicity amplitude approximatiorl^ ^ l ^^ l The numerical calculation 
beyond the helicity approximation differs only about a few percen|l24i 

By fitting the overall suppression for one nuclear target, we obtain 
the only parameter in our calculation, C{Q^) = 0.0060 GeV^ with 
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Fig. 26. Predicted nucl ear modification of jet fragmentation function is compared to 
the HERMES data ^HHI on ratios of hadron distributions between A and d targets in 
DIS. 
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Fig. 27. Energy dependence of the nuclear modification compared with the HERMES 
datalMI 



Q!s(Q^) — 0.33 at ~ 3 GeV^. This parameter is also related to nu- 
clear broade ning of the transverse momentum of the Drell-Yan dilept on in 
pA coUision^J^, (Ag^.) « Cnas/NcXA- With an experimental''^^ value 
of (Aql) = 0.016^1/3 GeV^ and a^iM^) = 0.21 (< M} >w 40 GeV^), 
one finds C{M?j) — 0.013 GeV^, which is about a factor 2 larger than the 
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value obtained in our fit to the HERMES data. The value of C determ ined 
from nuclear broadening in photo-production of a di-jet is even largeff^ 
at ~ 4py « 64 GeV^. Such a strong scale dependence of C{Q^) is in line 
with one's expe ctatio n, since it is related to gluon distribution xg{x, Q^) at 
small X in nuclePU. 

In principle the modification of the fragmentation functions would be the 
only experimental effect of induced gluon radiation via multiple scattering. 
One can never directly measure the energy loss of the leading quark. The 
net effect of the energy loss is the suppression of leading particles on one 
hand and the enhancement of soft particles on the other, leading to the 
modification of the fragmentation functions. One can then experimentally 
characterize the parton energy loss via the momentum transfer from large 
to small momentum regions of the fragmentation functions. 

Upon a close examination of Ea. Hll()|l . wc see that the first term is the 
renormalized fragmentation function in vacuum. The rest is particle produc- 
tion induced by the rescattering of the quark through the nuclear medium. 
In particular, the last term is particle production from the fragmentation 
of the gluon which is induced by the secondary scattering. Such particle 
production is at the expense of the energy loss of the leading quark. We 
can thus quantify the quark energy loss by the momentum fraction carried 
by the radiated gluon. 



{l\Zg){xB-,ti ) ^ j ~f~ j dz — zA-fq^gq{z,XB,XLjT) 



dz[l + (1 - ;z) 



(114) 



Using the approximation for the modified twist-four parton matrix elements 
in Ea. ljll^ll . we have 



{^zg){xB,^i')^c^^-^ dz 



Nc xaQ^ Jo z(l-z) 



X 



'^^?^(l_f)(l_e-iM), (115) 

^ 



where = ji^ /2p^q z(l — z) = xb/z{1 — z) if we choose the factorization 
scale as — Q^. When a;^ ^ a;_B <C 1 we can estimate the leading quark 
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energy loss roughly as 

{Azg){xB,fi )~C 50f[ln— -]. (116) 

^ X Xb ^ 

Since xa = lA/i?yi, the energy loss {^Zg) thus depends quadratically on 
the nuclear size. 

In the rest frame of the nucleus, = ttln, q" = j/, and xb = 
Q^/2p'^q~ = Q^/2mjvJ^. One can get the averaged total energy loss as 
AE = v{Azg) w q{Q'^)al{Q'^)mNR\{CAlN^)i\TL{l/2xB). With the de- 
termined value of C, {xb) ~ 0.124 in the HERMES experimenJinni and 
the average distance {La) — Ra\/2/'k for the assumed Gaussian nuclear 
distribution, one gets the quark energy loss dE/dL « 0.5 GeV/fm inside a 
Au nucleus. 



4.5. Energy Loss in Hot Medium at RHIC 

To extend our study of modified fragmentation functions to jets in heavy- 
ion collisions and to relate to results obtained in the opacity expansion 
approach, we can assume {k'^) « (the Debye screening mass) and a gluon 
density profile p{y) — {to/t)9{Ra ~ y)p() for a 1-dimensional expanding 
system. Since the initial jet production rate is independent of the final 
gluon density which can be related to the parton-gluon scattering cross 
sectior^l [asa;TG(xT) ^ A^^fg], one has then 

^fl^^^^p^l dyagp{y)[l-cos{y/rf)], (117) 

where Tf — 2Ez{\ — z^jl^, is the gluon formation time. One can recover 
the form of energy loss in a thin plasma obtained in the opacity expansion 
approach^fii 

{Azg)^^-^l\z f^du'^^'-^^' 



TT Jo Jo U{1 + U) 



/ dragpir) 

J To 



1 



2Ez{l - z) 



(118) 

Keeping only the dominant contribution and assuming Gg w Ca2iio?^l p? 
(Ca=l for qg and 9/4 for gg scattering), one obtains the averaged energy 
loss, 

,dE^ ^ i:CaGAO?s /-^"^ , 2E 

K-JTl- T. / «iTp(Tr- To) In — ^. (119) 

dL Ra Jro Tp^ 
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Neglecting the logarithmic dependence on t, the averaged energy loss in a 
l-dimensional expanding system can be expressed as 

dE ^dEo2ro 

where dEo/dL cx pqRa is the energy loss in a static medium with the same 
gluon density po as in a 1-d expanding system at time tq. Because of the 
expansion, the averaged energy loss {dE/dL)id is suppressed as compared 
to the static case and does not depend linearly on the system size. This 
could be one of the reasons why the effect of parton energy loss is found to 
be negligible in A A colhsions at = 17.3 Ge^CHl. 

An effective model of modified fragmentation functions was proposed in 
Ref.™ 

« Y^^a.. (y^) , (121) 

with Az to account for the fractional parton energy loss. This effective 
model is found to reproduce the pQCD result from Ec ljllUII very well, but 
only when Az is set to be Az w 0.6(zg). Therefore the actual averaged 
parton energy loss should be /S.E/E = 1.6Az with Az extracted from the 
effective model. The factor 1.6 is mainly caused by the unitarity correction 
effect in the pQCD calculation. A similar effect is also found in the opacity 

expansion approadS^l. 

Both PHENIX and STAR experiments have reportedly strong suppres- 
sion of high pt hadrons in central Au + Au collisions at -^i = 130 and 200 
GeV, indicating for the first time a large parton energy loss in heavy-ion 
collisions. To extract the parton energy loss, we compare the data with 
the calculated hadron pT spectra in heavy-ion collisions using the above 
effective model for medium modified jet fragmentation functions"^^. Shown 
in Fig. I2H1 are the nuclear modification factors Raa{pt) as the ratios of 
hadron spectra in AA {pA) and pp collisions normalized by the number 
of binary collisions^''. Parton shadowing and nuclear broadening of the in- 
trinsic kT are also taken into account in the calculation which decrsibes pA 
data for energies up to -Js = 40 Ge\^^. The nuclear fcr-broadening gives 
the Cronin enhancement at large pt in pA collisions, where there is no par- 
ton energy loss induced by a hot medium. Fitting the PHENIX data yields 
{dE/dL)id ~ 0.34(ln£;/ln5) GeV/fm, including the factor of 1.6 from the 
unitarity correction effect. We consider only tt^ data here, since at large pt 
the charged hadrons are dominated by baryons, which could be influenced 
mainly by non-perturbative dynamics 
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Fig. 28. Calculated nuclear modification factor of px spectra for d + Au and central 
Au + Au collisions at x/s = 130 (solid) and 200 GeV (dashed) as compared to PHENIX 
data El The lower dashed line uses the energy loss as given in Eg. 11221 . 



Taking into account the expansion, the averaged parton energy loss 
extracted from the PHENIX data would be equivalent to {dE/dL)o = 
0.34(i?yi/2ro) ln£'/ln5 in a static system with the same gluon density as 
the initial value of the expanding system at tq. With Ra ~ 6 fm and 
To ~ 0.2 fm, this would give {dE/dL)o w 7.3 GeV/fm for a 10-GeV parton, 
which is about 15 times of that in a cold Au nucleus, as extracted from 
the HERMES data. Since the parton energy loss is directly proportional 
to gluon density of the medium, this implies that the gluon density in the 
initial stage of Au + Au collisions at tq = 0.2 fm/c is about 15 times higher 
than that inside a cold nucleus. We can predict the 7r° spectra at ^/s = 200 
GeV as given by the dashed lines in Fig. |2H1 assuming that the initial par- 
ton density in central Au + Au collisions at ^/s = 200 GeV is about 10% 
higher than at 130 GeV. 

In Fig. EHl we have shown the predicted suppression factor at large pr 
with the same logarithmic energy-dependence of the energy loss as in the 
cold nuclear matter which gives a suppression factor that increases with px- 
However, new RHIC data'^' at 200 GeV show a constant suppression factor 
at larg er jP r- This might be the indication of the importance of the detailed 
balance^ which gives much stronger energy dependence of the energy loss. 
We show as the lower dashed line in Fig.EHlthe result for an effective parton 
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energy loss 

AS cx (£;/^- 1.6)^-27(7.5 + £;/Ai) (122) 

which is parameterized according to the result from where both 

stimulated gluon emission and thermal absorption are included in the cal- 
culation of the total energy loss. The detailed balance between emission and 
absorption both reduces the effective parton energy loss and increases the 
energy dependence. The threshold is the consequence of gluon absorption 
that competes with radiation that effectively shuts off the energy loss. The 
parameter fi is set to be 1 GeV in the calculation shown in Fig. EHI 

5. Summary 

In this report we reviewed two recent approaches to the problem of non- 
Abelian radiative energy energy loss in dense but finite QCD matter. In the 
first section, we highlighted some of the striking new high pT phenomena 
observed for the first time in A + A reactions at RHIC with ^/s > 100 
AGeV. An interesting pattern of hadron suppression, already beginning at 
moderate px > 5 GeV, was seen in the hadron flavor dependence of single 
inclusive spectra, in the large azimuthal asymmetry, and in the prelimi- 
nary two-hadron correlations. We interpret these phenomena as manifes- 
tations of jet quenching in ultra-dense matter produced in such reactions. 
Our predictions for these phenomena in both approaches are reviewed in 
later sections and depend on the energy loss, AE = J dx dE/dx, of fast 
quarks and gluons propagating through rapidly expanding QCD matter. 
The two approaches, GLV and WW/WOGZ reviewed here, provide a sys- 
tematic way to compute AE via an opacity or higher twist expansion in 
finite nuclear ma tter. El sewhere reviewed asymptotic approaches such as 
the BDMS/Z/sW^Smni 

are designed for applications to "thick" or macro- 
scopic media at asymptotic energies. An analytic approximation to the sum 
of the infinite multiple collision series is obtained through an approximate 
color dipole quantum diffusion analogous to the Moliere series in electro- 
dynamics. The c omplications due to finite kinematic bounds are neglected. 
As shown i the phenomenological applications of the asymptotic ex- 
pressions tend to overpredict quenching at RHIC and lead to a too rapid 
variation of the suppression factor with px, inconsistent with the RHIC 
data. In our approach, on the other hand, the opacity series is computed to 
arbitrary order in opacity for applications to finite opacity systems where 
the non-Gaussian (Rutherford) tails of distributions are not yet eclipsed by 
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the approximate Gaussian small pt component. In addition, our expres- 
sions can be applied to arbitrary 3D expanding and time dependent media 
such as created in nuclear collisions. 

The GLV reaction operation approach discussed in Section 2 is based 
on a general algebraic recursive method that describes the propagation and 
interaction of systems through dense nuclear matter, taking into account the 
severe destructive LPM interference effects due to the long formation times 
of ultra-relativistic partons. It provides a means to compute the multiple 
collision amplitudes to any order x" in the opacity x — J dz ap. So far it has 
been successfully applied to obtain solutions for th e nuclear broadeninj^ 
and the final state medium induced radiatiorP^lSl resulting from multiple 
elastic and inelastic projectile scatterings. The range of applicability of the 
calculations can be significantly extended by careful treatment of kinematic 
bounds. 

The final state double differential distributions of jets and gluons are 
presented as an infinite series in powers of the mean number of scatter- 
ings X- For elastic scatterings this series can be summed to all orders (or 
equivalently all twist parton-parton correlations) to reproduce the standard 
Glauber theory result, but more general than the simple color dipole ap- 
proximation used in asymptotic analyses. For inelastic processes, a closed 
sum to all orders can unfortunat ely on ly be carried out in the dipole approx- 
imation as in BDMS/Z/SW (se a '^ 1 for most recent developments). How- 
ever, our analyti c ^^ '^^' e xpressions at any finite order in opacity can be eval- 
uated numericalljl^^^^. Numerical evaluation of the expressions through 
the first three orders (up to twist 8) has been carried out for phenomeno- 
logical applications. Each power in opacity adds a twist 2 parton-parton 
correlation. Future work via this approach includes computation of multi- 
ple gluon emission beyond the Poisson approximation and the broadening 
and radiation of dipole-like, possibly heavy, qq systems. Applications to 
heavy quark energy loss including both the Ter-Mikayelian gluon disper- 
sion effects as well as the "dead cone" effect are also underway . 

The WW approach reviewed in Section 3 extends the calculation of 
energy loss to include the positive feedback (jet acceleration) due to ab- 
sorption of thermal gluons in the medium. Absorption counteracts the in- 
duced energy loss for jet momenta less than the typical thermal energy scale 
^ ST. In an expanding hydrodynamic medium with transverse boost rapid- 
ity rjT we expect that this absorption feedback contribution is blue shifted 
to higher momenta pr ~ 3Te''^. Thus jet quenching cannot suppress the 
spectrum below the local equilibrium hydrodynamic limit. 
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In Section 4, several recent works were combined (referred to as the 
WOGZ approach*^^ '^^ 104^ elaborate on a general twist-expansion 
of multiple parton scattering beyond the GW model of the medium. In 
that approach one can calculate explicitly the modified parton fragmenta- 
tion function up to twist 4 thus far. The LPM interference effect is then 
embedded in the twist-four parton matrix elements of the nucleus which 
also contains the fundamental properties of the nuclear medium - parton 
density correlations inside a nucleus. These matrix elements replace the De- 
bye screened interactions used in the GLV approach. One can demonstrate 
explicitly that the quadratic dependence of the modification of fragmen- 
tation functions and the effective parton energy loss on the nuclear size 
Ra is caused both by the LPM interference and the specific form of gluon 
radiation spectra in QCD. The predicted nuclear modification of the frag- 
mentation function, both the energy and nu clear dependence, is found to 
agree well with the recent experimental dataPSI 

on jet fragmentation in 

e + A. This is an important test of both the twist or opacity expansions 
since it shows the dominance of the twist 4 (first order in opacity) contri- 
bution to the energy loss in finite nuclear systems as found numerically in 
the GLV expansion up to twist 8. Extending the results to a parton prop- 
agating in a hot QCD medium, we have shown that the twist expansion 
to order 4 is equivalent to the first order opacity GLV result under certain 
simplifying assumptions about the form of the twist 4 matrix elements. The 
phenomenological apphcation of WOGZ to both DIS on cold nuclear targets 
and high-energy heavy-ion collisions suggests that the parton energy loss 
in an expanding system at RHIC would be equivalent to {dE/dx)o ~ 7.3 
GeV/fm in a static medium, which is almost 15 times higher than that in 
a cold Au nucleus. 

If the jet quenching pattern is confirmed by further measurements and 
theoretical refinements, current RHIC data may have already provided the 
first tomographic evidence that initial parton densities on the order of 100 
times nuclear matter density were achieved in Au + Au collisions. The full 
analysis of the flavor composition, shape, and azimuthal moments of the 
high pt spectra appears to be a promising diagnostic probe of the evolution 
of the produced quark-gluon plasma. 

However, in spite of the consistency of our jet tomography analysis with 
current RHIC data, it is still too early to draw definitiv e con clusions. The 
main uncertainty is the magnitude of gluon shadowin^^^l in the initial 
nuclear wavefunction. While estimates of shadowing vary wildly, unfortu- 
nately nothing is yet known experimentally on this important question. 
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QCD analysijni suggests that non-Unear corrections to DGLAP evolution 
may cause the gluon densi ty to saturate at small x. Saturation is found in 
classical Yang-Mills j^Q^jglJIininnillZinHl as weH. in Ref.^^ it was pro- 
posed that gluon saturation or shadowing phenomena alone may account 
for a significant part of the observed high pT hadron suppression pattern. 

AO 

Our estimates of shadowing bas ed on ES K'9&^ do not sup port this idea 
and neither do recent calculation^^^^^in the MQ approach !^ of the up- 
per limit of the saturation scale and the rate of disappearance of its effects 
on the small and moderate pt spectra. Fortunately, a decisive experimental 
test is now under way at RHIC via d+Au reactions. If our theory is correct, 
and the observed suppression pattern in Au + Au is due to jet quenching 
via final state interactions in the dense QCD matter produced, then instead 
of suppression, the Cronin effect is predicted to dominate over shadowing 
in the x > 0.01 range accessible at RHIC. We therefore predict as an up- 
per limit an enhancement of the moderate pt hadrons by ^ 10 — 30% in 
d + A relative to binary scaled p + p. On the other hand, if the saturation 
model is correct, and the suppression pattern in Au -\- Au is a consequence 
of deep gluon shadowing in the initial state, then a suppression a bout 30% 
of moderate and high px hadrons is predicted to occur in ci -I- 

The experimental answer to this decisive question may be known by the 
publication time of this review. In either case, the answer will be exciting. 

Acknowledgments 

M.G. and I.V. gratefully acknowledge extensive collaboration with Pe- 
ter Levai on the GLV approach. X.N.W. and B.W.Z. would like to thank 
Enke Wang on extensive collaboration on the WW approach with detailed 
balance and many fruitful discussions. X.N.W. would also like to acknowl- 
edge collaborations with X. F. Guo and J. Osborne on the twist expansion 
approach. This work was supported by the Director, Office of Energy Re- 
search, Office of High Energy and Nuclear Physics, Division of Nuclear 
Physics, and by the Office of Basic Energy Science, Division of Nuclear 
Science, of the U.S. Department of Energy under Contracts No. DE-FG02- 
87ER40371, DE-FG02-93ER40764, and DE-AC03-76SF00098 and by Na- 
tional Natural Science Foundation of China under project No. 19928511 
and No. 10135030. 

References 

1. K. Adcox et al. [PHENIX Collaboration] , |axXiv:nucl-ex/0207009| 



February 6, 2008 21:23 WSPC/Trim Size: 9m x 6in for Review Volume GVWZ 



66 



2. C. Adler et ai, Phys. Rev. Lett. 89, 202301 (2002). 

3. S. Mioduszewski [PHENIX Collaboration], arXiv:nucl-ex/021002l] 

4. D. d'Enterria [PHENIX Collaboration], arXiv:hep-ex/020905lJ 

5. G. J. Kunde [STAR Collaboration], arXiv:nucl-ex/0211018. 

6. C. Adler et al. [STAR Collaboration], arXiv:nucl-ex/0210'^ 

7. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 907032301 (2003). 

8. K. Filimonov [STAR Collaboration], arXiv:nucl-ex/0210027 

9. P. Levai, G. Papp, G. Fai, M. Gyulassy, G. G. Barnafoldi, L Vitev and 
Y. Zhang, Nucl. Phys. A 698, 631 (2002). 

10. Proc. 15th Int. Uhra. Nucleus-Nucleus Coll. (Quark Matter 2001, Jan. 
SUNY), Nucl.Phys. A 698 (2002) Ic; Proc. 16th (Quark Matter 2002, July 
Nantes, France); Nucl. Phys. A in press. 

11. J. D. Bjorken, FERMILAB-PUB-82-59-THY and erratum (unpublished) 

12. M. H. Thoma and M. Gyulassy, Nucl. Phys. B 351, 491 (1991). 

13. M. H. Thoma, Phys. Rev. D 49, 451 (1994). 

14. D. A. Appel, Phys. Rev. D 33, 717 (1986). 

15. J. P. Blaizot and L. D. McLerran, Phys. Rev. D 34, 2739 (1986). 

16. M. Rammerstorfer and U. W. Heinz, Phys. Rev. D 41, 306 (1990). 

17. M. Gyulassy and M. Plumer, Phys. Lett. B 243, 432 (1990). 

18. M. Gyulassy and M. Plumer, Nucl. Phys. A 527, 641 (1991). 

19. M. Gyulassy, M. Plumer, M. Thoma and X. N. Wang, Nucl. Phys. A 538, 
37C (1992). 

20. X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 (1991). 

21. X. N. Wang and M. Gyulassy, Phys. Rev. D 45, 844 (1992). 

22. M. Gyulassy and X. N. Wang, Comput. Phys. Commun. 83, 307 (1994). 

23. X. N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 (1992). 

24. M. M. Aggarwal et al. [WA98 Collaboration], Phys. Rev. Lett. 81, 4087 
(1998) [Erratum-ibid. 84, 578 (1998)]. 

25. X. Wang and M. Gyulassy, Phys. Rev. Lett. 86, 3496 (2001). 

26. E. Wang and X. N. Wang, Phys. Rev. C 64, 034901 (2001) 

27. C. Adler et at [STAR Collaboration], Nucl. Phys. A 698, 64 (2002). 

28. W. A. Zajc et at [PHENIX Collaboration], Nucl. Phys. A 698, 39 (2002). 

29. G. Roland [PHOBOS Collaboration], Nucl. Phys. A 698 54 (2002). 

30. I. G. Bearden et al. [BRAHMS Collaboration], Nucl. Phys. A 698, 29 (2002). 

31. P. F. Kolb, P. Huovinon, U. W. Heinz and H. Heiselberg, Phys. Lett. B 500, 
232 (2001). 

32. P. Huovinen, P. F. Kolb, U. W. Heinz, P. V. Ruuskanen and S. A. Voloshin, 
Phys. Lett. B 503, 58 (2001). 

33. P. F. Kolb, U. W. Heinz, P. Huovinen, K. J. Eskola and K. Tuominen, Nucl. 
Phys. A 696, 197 (2001). 

34. D. Teaney, J. Lauret and E. V. Shuryak, [iiucl-th /01040'4T1 

35. X. N. Wang, Phys. Rev. C 63, 054902 (20017 

36. M. Gyulassy, I. Vitev and X. N. Wang, Phys. Rev. Lett. 86, 2537 (2001). 

37. M. Pluemer, M. Gyulassy and X. N. Wang, Nucl. Phys. A 590, 511C (1995). 

38. J. F. Owens, Rev. Mod. Phys. 59 (1987) 465. 

39. X. N. Wang, Phys. Rev. C 61, 064910 (2000). 



February 6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



JET QUENCHING AND RADIATIVE ENERGY LOSSIN DENSE NUCLEAR MATTERS! 

40. I. Vitev, arXiv:nucl-th/0302002 

41. K. J. Eskola and H. Honkanen, Nucl. Phys. A 713 167 (2003). 

42. I. Vitev, arXiv:hep-ph/0212109 

43. X. N. Wang, Phys. Rev. C 58, 2321 (1998). 

44. X. N. Wang, Z. Huang and I. Sarcevic, Phys. Rev. Lett. 77, 231 (1996). 

45. X. N. Wang and Z. Huang, Phys. Rev. C 55, 3047 (1997). 

46. R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, JHEP 0109, 033 
(2001). 

47. M. Gyulassy, P. Levai and I. Vitev, Phys. Lett. B 538, 282 (2002). 

48. K. J. Eskola, V. J. Kolhinen and C. A. Salgado, Eur. Phys. J. C 9, 61 (1999). 

49. L. D. Landau and I. Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 92, 535 
(1953). 

50. L. D. Landau and I. Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 92, 735 
(1953). 

51. A. B. Migdal, Phys. Rev. 103, 1811 (1956). 

52. M. Gyulassy, P. Levai and I. Vitev, Nucl. Phys. B 594, 371 (2001). 

53. M. Gyulassy, P. Levai and I. Vitev, Phys. Rev. Lett. 85, 5535 (2000). 

54. M. Gyulassy, P. Levai and I. Vitev, Nucl. Phys. B 571, 197 (2000). 

55. M. Gyulassy, P. Levai and I. Vitev, Nucl. Phys. A 661, 637 (1999). 

56. R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. 
Phys. B 483, 291 (1997). 

57. R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. 
Phys. B 484, 265 (1997) 

58. R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, Nucl. Phys. B 531, 
403 (1998) 

59. R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, Phys. Rev. C 60, 
064902 (1999) 

60. B. G. Zakharov, JETP Lett. 63, 952 (1996). 

61. B. G. Zakharov, JETP Lett. 65, 615 (1997). 

62. B. G. Zakharov, JETP Lett. 73, 49 (2001). 

63. B. G. Zakharov, JETP Lett. 76, 201 (2002). 

64. U. A. Wiedemann, Nucl. Phys. A 690, 731 (2001). 

65. U. A. Wiedemann, Nucl. Phys. B 588, 303 (2000). 

66. U. A. Wiedemann, Nucl. Phys. B 582, 409 (2000). 

67. C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 89, 092303 (2002) 

68. R. Baier, D. Schiff and B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50, 37 
(2000). 

69. X. N. Wang and X. F. Guo, Nucl. Phys. A 696, 788 (2001). 

70. X. F. Guo and X. N. Wang, Phys. Rev. Lett. 85, 3591 (2000). 

71. E. Wang and X. N. Wang, Phys. Rev. Lett. 89, 162301 (2002). 

72. M. Gyulassy and X. Wang, Nucl. Phys. B 420, 583 (1994). ; Wang X N, 
Gyulassy M, Plumer M, Phys. Rev. D 51 (1995) 3436. 

73. M. Gyulassy, I. Vitev, X. N. Wang and P. Huovinen, Phys. Lett. B 526, 
301 (2002). 

74. M. Gyulassy, P. Levai and 1. Vitev, Phys. Rev. D 66, 014005 (2002). I. Vitev 
and J.W. Qiu, work in preparation. 



6,2008 21:23 WSPC/Trim Size: 9m x 6in for Review Volume GVWZ 



68 



75. I. Vitev and M. Gyulassy, Phys. Rev. Lett. 89, 252301 (2002). 

76. A. Accardi, arXiv:hep-ph/0212148 

77. Y. Zhang, G. Fai, G. Papp, G. G. Barnafoldi and P. Levai, Phys. Rev. C 
65, 034903 (2002). 

78. B. Z. Kopehovich, J. Nemchik, A. Schafer and A. V. Tarasov, Phys. Rev. 
Lett. 88, 232303 (2002). 

79. F. Arleo, Phys. Lett. B 532, 231 (2002). 

80. J. Velkovska [PHENIX collaboration], Nucl. Phys. A 698, 507 (2002). 

81. K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 88, 242301 
(2002). 

82. N. Xu and M. Kaneta, Nucl. Phys. A 698, 306 (2002). 

83. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 87, 262302 (2001). 

84. K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 89, 092302 
(2002). 

85. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 89, 092301 (2002) 
arXiv:nucl-ex/0203016 . 

86. 1. Vitev and M. Gyulassy, Phys. Rev. C 65, 0419 02 (2002). 

87. 1. Vitev and M. Gyulassy, |arXiv:hep-ph/0208108| 

88. 1. Vitev, M. Gyulassy and P. Levai, 'arXiv:hep-ph/0109198" 

89. G. C. Rossi and G. Veneziano, Nucl. Phys. B 123, 507 (1977). 

90. L. Montanet, G. C. Rossi and G. Veneziano, Phys. Rept. 63, 149 (1980). 

91. D. Kharzeev, Phys. Lett. B 378, 238 (1996). 

92. S. E. Vance and M. Gyulassy, Phys. Rev. Lett. 83, 1735 (1999). 

93. S. E. Vance, M. Gyulassy and X. N. Wang, Phys. Lett. B 443, 45 (1998). 

94. J. Y. OUitrauh, Phys. Rev. D 46, 229 (1992). 

95. E. V. Shuryak, Phys. Rev. C 66, 027902 (2002). 

96. T. Hirano and Y. Nara, arXiv:nucl-th/0301042 

97. T. Hirano and Y. Nara, arXiv:nucl-th/0211096 

98. E. Wang and X. N. Wang, Phys. Rev. Lett. 87, 142301 (2001). 

99. E. Braaten and R. D. Pisarski, Nucl. Phys. B 337, 569 (1990). 

100. V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); Yu. 
L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977); G. Altarelh and G. Parisi, 
Nucl. Phys. B126, 298 (1977); 

101. J. Osborne and X.-N. Wang, Nucl. Phys. A 710, 281 (2002) 
arXiv:hep-ph/0204046 . 

102. M. Luo, J. Qiu and G. Sterman, Phys. Lett. B279, 377 (1992); M. Luo, J. 
Qiu and G. Sterman, Phys. Rev. D50, 1951 (1994); M. Luo, J. Qiu and G. 
Sterman, Phys. Rev. D49, 4493 (1994). 

103. J. Qiu, Phys. Rev. D 42, 30 (1990). 

104. B. W. Zhang and X.-N. Wang, arXiv:hep-ph/0301195 

105. A. H. Mueller and J. Qiu, Nucl. Phys. B268, 427 (1986). 

106. A. Airapetian et al. [HERMES Collaboration], Eur. Phys. J. C 20, 479 
(2001); V. Muccifora [HERMES Collaboration], arXiv:hep-ex/0106088 

107. X. F. Guo, Phys. Rev. D 58, 114033 (1998). 

108. P. L. McGaughey, J. M. Moss and J. C. Peng, Ann. Rev. Nucl. Part. Sci. 
49, 217 (1999); J. C. Peng, priviate communication. 



6, 2008 21:23 



WSPC/Trim Size: 9in x 6in for Review Volume 



GVWZ 



JET QUENCHING AND RADIATIVE ENERGY LOSSIN DENSE NUCLEAR MATTERm 

109. X.-N. Wang, Phys. Rev. Lett. 81, 2655 (1998). 

110. Note Added in Proof: On completion of this report an instructive new analy- 
sis comparing the BDMS/Z/SW and the GLV approaches became available: 
C. A. Salgado and U. A. Wiedemann, arXiv:hep-ph/0302184,, U. A. Wiede- 
mann, arXiv:hep-ph/0302194 

111. F. Arleo, JHEP 0211, 044 (2002). 

112. M. Djordjevic and M. Gyulassy, arXiv:nucl-th/0302069| 

113. N. Armesto and C. A. Salgado, arXiv:hep-ph/0301200 

114. A. H. Mueller and J. w. Qiu, Nucl. Phys. B 268, 427 (1986). 

115. L. D. McLerran and R. Venugopalan, Phys. Rev. D 50, 2225 (1994). 

116. L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 3352 (1994). 

117. E. lancu, A. Leonidov and L. McLerran, arXiv:hep-ph/0202270 

118. A. H. Mueller, iarXiv:hep-ph/011124'4l 

119. D. Kharzeev, E. Levin and L. McLerran, "arXiv:hep-ph/0210332] 

120. K. J. Eskola, H. Honkanen, V. J. Kolhinen, J. W. Qiu and C. A. Salgado, 
arXiv:hep-ph/0211239 

121. K. J. Eskola, H. Honkanen, V. J. Kolhinen, J.W. Qiu and C. A. Salgado, 
|arXiv:hep-ph/03 02185 



